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GENERAL  THEORY  OF  INTERMOLECULAR  FORCES 


A.D.  Buckingham 
University  Chemical  Laboratory 
Lensfield  Road 
Cambridge 
CB2  1EW 
England 


ABSTRACT 


There  is  much  current  interest  both  in  inicr-aml  intra-inolecuiar  potential  energy 
surfaces.  The  structure  and  properties  of  van  der  Waals  molecules  and  clusters  provide 
an  important  source  of  information  about  molecular  interactions.  This  information  can 
be  used  to  generate  intermolecular  potentials  that  can  be  useful  in  descriptions  of  larger 
systems,  such  as  condensed  phases,  solutions  of  surfactants,  and  biomacromolecules. 
It  is  convenient  to  divide  intermolecular  potentials  into  long-range  and  short-range 
components.  The  former  are  related  via  perturbation  theory  to  the  charge  distribution 
and  polarizabilities  of  the  free  molecules,  and  the  resulting  long-range  potentials  vary  as 
an  inverse  power  of  the  separation  between  the  molecules.  The  short-range  interactions 
result  from  the  overlap  of  the  electron  clouds  of  the  interacting  molecules  and  diminish 
exponentially  with  the  separation.  By  dividing  the  molecules  into  atoms  or  small 
groups  of  atoms,  it  is  possible  to  obtain  convenient  and  convergent  representations  of 
the  potential.  The  approach  can  be  used  to  provide  a  theoretical  basis  for  the  popular 
site-site  potential  models  that  are  now  used  extensively  in  computer  simulations. 
Attention  will  be  paid  to  the  additivity  or  non-additivity  of  potentials,  to  the  role  of  the 
solvent  in  solutions,  and  to  changes  in  the  electronic  properties,  such  as  the  dipole 
moment,  that  result  from  intermolecular  forces. 
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A  VIBRATIONAL  MOLECULAR  FORCE  FIELD 
FOR  MACROMOLECULAR  MODELLING 


Gerard  VERGOTEN 
FacuitedePharmacie,  INSERM  U279 
Universite  de  LILLE  II  -  LILLE  -  FRANCE 


An  intramolecular  force  field  has  been  developped  for  various  building  blocks  with 
biological  interest!  amino-acids,  peptides,  nucleotides  and  saccharides  ) .  A  UREY-  BRAD  LEY’ - 
SHIMANOUCHI  type  expression  has  been  used  for  the  potential  energy  function.  The 
correspondingparametersareobtaiaedontfaebasis  of  vibrational  frequencies  through  normal  modes 
analysis  for  the  crystals.  In  order  to  extract  the  best  sec  among  the  various  force  fields  able  to 
reproduce  frequencies,  additional  data  are  needed.  At  the  present  time  they  consist  of  infrared  and 
resonance  Ramanintensities  and  temperature  factors  given  by  the  X  rays  studies.  The  transferability 
of  the  final  set  is  also  checked  on  larger  systems  and  then  used  for  molecular  mechanics  and 
dynamics  calculations.  A  harmonic  dynamics  study  of  bovine  pancreatic  trypsin  inhibitor  and  a 
turn  of  the  left  handed  Z-  DNA  heiix  will  be  presented  and  compared  with  previous  published  works. 
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MODELLING  OF  BIOMOLECULAR  PROPERTIES  BY  MOLECULAR  DYNAMICS  SIMULATION 


W.F.  van  Gunsteren 

Department  of  Physical  Chemistry 
University  of  Groningen 
Nijenborgh  16 
4  7  AG  Groningen 
The  Netherlands 

and 

Department  of  Physics 
Free  University 
P.0.  Box  7161 
1007  MC  .Amsterdam 
The  Netherlands 


During  the  last  decade  it  has  become  feasible  to  simulate  the  dynamics 
of  biologically  relevant  molecules  on  a  computer.  The  method  of  molecular 

dynamics  (MD)  solves  Newton’s  equations  of  motion  for  a  molecular  system, 
which  results  in  trajectories  for  all  atoms  in  the  system.  From  these 
atomic  trajectories  a  variety  of  properties  can  be  calculated. 

Computer  simulation  of  molecular  systems  aims  at  macroscopic  behaviour 
to  be  computed  from  microscopic  interactions.  The  main  contributions  of  the 

microscopic  point  of  view  are  (1)  understanding,  (2)  interpretation  of 
experimental  results,  (3)  semiquantitative  estimates  of  experimental 
results  and  (4)  the  capability  to  interpolate  or  extrapolate  experimental 
data  into  regions  that  are  hard  to  reach  in  the  laboratory. 

One  of  the  two  basic  problems  in  the  field  of  molecular  modelling  and 
simulation  is  hosv  to  efficiently  search  the  vast  configuration  space  which 
is  spanned  by  all  possible  molecular  conformations  for  the  global  low 
(free)  energy  regions  which  will  be  populated  by  a  molecular  system  in 
thermal  equilibrium.  The  other  basic  problem  is  the  derivation  of  a 
sufficiently  accurate  interaction  energy  function  or  force  field  for  the 

molecular  system  of  interest.  It  belongs  to  the  art  of  computer  simulation 
to  choose  the  unavoidable  assumptions,  approximations  and  simplifications 
of  the  molecular  model  and  computational  procedure  such  that  their 
contributions  to  the  overall  inaccuracy  are  of  comparable  size,  without 
affecting  significantly  the  property  of  interest. 

Some  practical  applications  of  computer  simulation  in  the  field  of 

(bio)chemistry  will  be  reviewed. 
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THEORETICAL  STUDY  OF  A  CONFORMATIONAL  CHANGE  IN  AN  ENZYME  : 


METHODOLOGY  AND  FIRST  RESULTS  ON  CITRATE  SYNTHASE 
Jean  Dump 

Laboratoire  de  Physique  Quantique  (I'A  505  du  C.N.R.S.) 
118,  route  de  Narbonne,  31062  TOULOUSE  CEDEX 


Following  the  lines  and  using  the  codes  developped  in  Harvard  by  M.  Karplus 
and  his  colleagues,  we  introduced  an  alternative  strategy  for  dynamical  calculations 
on  proteins  at  various  time  scales.  It  essentially  consists  of  the  following  steps. 

(i)  We  construct  a  set  of  internal  vectors  bv  the  principle  of  the  binary  tree, 
adapted  to  the  peculiar  structure  of  each  residue  and  to  the  secondary  and  tertiary 
structures  of  the  protein.  The  transformat  ion  from  Cartesian  coordinates,  velocities 
and  energy  gradients  to  the  new  set  and  vice  versa  needs  purely  topological,  cons¬ 
tant  matrices  ;  the  kinetic  energy  remains  diagonal  and  thus  simple  Newtonian  dyna¬ 
mics  may  be  used. 


(ii)  From  these  internal  vectors  we  generate  a  set  of  relative  polar  coordi¬ 
nates  and  determine  through  a  dynamical  test  the  frequencies  associated  with,  each 
of  these  coordinates,  which  permits  a  classification  of  slow  and  fast  modes. 

(iii)  We  then  compared  a  regular  dynamics  by  the  CHARMM  programs,  using  a 

1  fs  integration  step,  with  a  constrained  dynamics  using  10  fs  steps  after  freezing 
all  modes  of  frequencies  larger  than  'O''  s  '.  It  appears  that  the  main  features  of 
the  dynamics  of  the  slower  modes  are  preserved,  apart  from  phase  shifts  occurinq 
during  the  freezing  process. 


( iv ;  Iteration  of  this  procedure  to 
averaged  forces  ire  used  for  the  hydrngc 


arger  time  scales  is  possible  prov 
bonds.  One  may  then  hope  t'  obtain 


ficant  collective  modes  by  diagonalizacion  of  the  correlation  matrix  of  slow  mode 
amp  1 i tudes . 
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Design  and  Synthesis  of  a  new  bifunctional  microprotein 
by  interactive  modelling 

D.  Le  Nguyen  A.  Heitz.  L.  Chiche.  D.  Nalis  and  B.  Castro 
Centre  CNRS-INSERM  de  Pharmaco-Endocrinologie 
34094  Montpellier.  France 

A  Fa'/el.  H  Mattras  and  M  A,  Coletti-Previero 
INSERM  U58.  rue  de  Navacelle. 

34000  Montpellier.  France 

C.  Gaboriau  and  J.P.  Mornon 

Laboratoire  de  Cristallographie.  Universite  Pierre  et  Marie  Cure. 

Place  Jussieu.  75005  Paris.  France 

EETI  II.  a  28  residues,  three  disulfide  bridges  microprotein  has  been  isolated  from  the 
cucurbitacaea  Ecballium  elaterium.  It  behaves  as  a  powerful  trypsin  inhibitor  with  a 
dissociation  constant  of  10I2M'.  Its  chemical  synthesis  has  been  achieved  that 
provided  enough  material  to  make  a  complete  2D  'H  NMR  study  :  then  the  3D  model 
was  built  using  the  DISGEO  method,  giving  the  tertiary  structure  and  the  disulfide 
bridges  arrangement.  The  DISGEO  structure  was  then  optimized  by  Restrained 
Molecular  Dynamics.  A  radiocrystallographic  study  of  the  EETI  11/  Trypsin  complex 
confirmed  the  NMR  attributions.  The  synthesis  of  peptides  bearing  selected  modifica¬ 
tions  at  the  reactive  site  has  been  successful.  The  new  products  show  high  inhibiting 
power  towards  elastase  or  chymotrypsin.  while  their  affinity  towards  trypsin  vanished. 
The  topological  analogy  of  EETI  II  to  the  known  Carboxypeptidase  A  inhibitor  from 
potato  suggested  to  build  a  chimeric  molecule  bearing  two  active  sites  against  both 
trypsin  and  carboxypeptidase  A.  A  32-peptide  designed  by  interactive  model  building 
was  synthetized  that  could  be  proven  to  bind  both  enzymes. 


t 


A 


J. 
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POTENTIEL  MAPS  OF  METHANE,  WATER  AND  METHANOL  IN  SILICALITE 
F.  VIGNE-MAEDER 

Institut  de  Recherches  sur  la  Catalyse 
2  avenue  Albert  Einstein 
69626  VILLEURBANNE  CEDEX 

Silicalite  is  a  form  of  pure  crystalline  silica  that 
represents  the  end  member  (Si/Aloo)  of  the  high-silica  ZSM-5 
zeolites.  It  provides  an  example  of  a  microporous  adsorbent 
without  any  adsorption  sites  of  chemical  nature,  which  adsorption 
properties  are  only  related  to  the  channel  structure. 

The  interaction  energy  of  silicalite  clusters  with  methane, 
methanol  and  water  has  been  calculated  and  optimized  with  respect 
to  the  orientation  for  many  positions  inside  the  channels.  The 
calculations  were  performed  in  the  atom-atom  approximation  with 
parameters  deduced  from  ab-initio  studies  of  small  systems.  The 
electrostatic  terms  are  evaluated  by  using  multicentered 
multipole  expansions,  limited  to  charges  and  dipoles  distributed 
on  all  the  atoms.  For  the  silicalite  crystal,  the  whole  multipole 
expansion  was  reconstructed  from  those  of  fragments,  especially 
dimers  (HO^SiOSitOH)^  and  monomers  Si(0H}^,  by  superposition  and 
substraction  of  parts  in  excess.  Concerning  the  silicalite 
clusters,  we  selected  a  set  of  cylinders  surrounding  the 
different  parts  of  the  channel  network.  All  silica  tetrahedra 
located  inside  one  of  these  cylinders  were  taken  into  account 
and  the  free  valences  were  saturated  with  hydrogen  atoms.  The 
potential  maps  represent  the  distribution  of  the  interaction 
energy  of  one  adsorbed  molecule,  averaged  over  different 

orientations . 

Approximate  values  of  diffusion  barriers  are  deduced  from 
the  potential  maps  and  the  initial  heats  of  adsorption  are 
evaluated  by  average  over  the  positions  and  the  orientations, 
with  the  hypothesis  of  a  3oltzman  distribution. 

Our  values  are  qualitatively  in  accordance  with  the 
experimental  results.  In  particular  the  orders  CH^  <  H^C  <  CH^OH 
for  the  heat  of  adsorption  and  CH^  >  H^O  >  CH^OH  for  the  mobility 
inside  the  crystal  are  well  reproduced.  The  values  of  the  initial 
heat  of  adsorption  (9-4  kcal/mol  for  CH^ ,  12.5  for  H.,0, 

16.7  for  CH^OH)  are  too  large  by  2  kcal/mol.  They  are  very 
sensitive  to  the  experimental  geometry  adopted  for  silicalite. 


7 


Chemisorption  on  Transition  Metal  Surfaces  and  Heterogeneous 


talysis:  Back-of-the-Envelope  Theoretical  Modeling 


Evgeny  Shustorovich 
Corporate  Research  Laboratories 
Eastman  Kodak  Company,  Rochester,  N.Y.  14650-2001,  U.S-A. 


The  lecture  will  concentrate  on  realistic  calculations  of  reaction  energy  profiles  on  metal 
surfaces,  comprising  heats  of  adsorbate  chemisorption  Q  and  activation  barriers  AE*  for 
adsorbate  transformations.  The  framework  is  the  ”back-of-the-envelope"  BOC-MP  (bond- 
order  conservation  Morse-potential)  model  providing  the  rigorous  analytic  formalism  to 
calculate  the  values  of  Q  and  AE*  for  diatomic  admolecules  (1).  This  formalism  can  be 
applied  to  polyatomic  molecules  as  well  if  their  total  bond  energies  in  the  gas-phase  are 
effectively  partitioned  into  two-center  contributions.  The  partitioning  procedure  resem¬ 
bling  the  use  of  Pauling  resonance  structures  will  be  described.  As  examples,  the  BOC-MP 
projections  concerning  CO  hydrogenation,  C2  hydrocarbon  transformations,  and  HCOOH 
decomposition  on  transition  metal  surfaces  will  be  discussed  in  detail,  in  particular  the 
preferred  coordination  modes  of  polyatomic  adsorbates  and  the  relative  preference  of 
various  conceivable  reaction  pathways.  We  will  demonstrate  the  basic  consistency  of  the 
BOC-MP  results  with  experiment,  especially  when  they  are  at  variance  with  other  the¬ 
oretical  approaches.  It  appears  that  the  BOC-MP  modeling  provides  a  practical  means  for 
understanding  and  projecting  surface  reactivity. 


(1)  E.  Shustorovich,  Surf.Sci.Repc.,  6(l9S6)  1-63;  Acc.Chem.Res.,  21(1988)  183-189. 
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Theoretical  models  for  intermolecular  potentials 
Anthony  J.  Stone 

University  Chemical  Laboratory,  Lensfield  Road,  Cambridge.  England 

It  is  now  possible  to  carry  out  accurate  ab  initio  calculations  on  molecular  com¬ 
plexes  by  a  variety  of  techniques.  The  supermolecule  approach  is  widely  used,  and  is 
capable  of  high  absolute  accuracy,  but  it  is  subject  to  Basis  Set  Superposition  Error, 
especially  when  electron  correlation  is  taken  into  account,  and  this  is  a  difficulty  when 
accurate  calculations  of  small  interaction  energies  are  required.  Perturbation  theory 
gives  a  more  detailed  description  of  the  interaction  than  the  supermolecule  approach, 
and  consequently  provides  more  physical  insight  into  the  nature  of  the  interaction,  but 
it  does  not  converge  if  the  interaction  is  too  strong,  that  is  when  the  distance  of  ap¬ 
proach  is  too  short,  although  in  practice  the  results  are  satisfactory  to  distances  inside 
the  potential  minimum.  Both  of  these  methods  require  calculations  to  be  carried  out 
at  a  wide  range  of  dimer  geometries  if  a  full  description  of  the  potential  energy  surface 
is  needed,  and  this  is  extremely  time-consuming. 

A  useful  alternative  approach  is  to  isolate  the  components  of  the  perturbation 
expansion,  namely  the  repulsion,  the  electrostatic  interaction,  the  induction,  the  dis¬ 
persion  and  the  charge  transfer  terms,  and  to  calculate  each  of  them  independently 
by  the  most  appropriate  technique.  Thus  the  electrostatic  interaction  can  be  calcu¬ 
lated  accurately  from  distributed  multipole  descriptions  of  the  individual  molecules, 
while  the  induction  and  dispersion  contributions  may  be  derived  from  polarizabilities 
of  the  individual  molecules.  This  approach  has  the  advantage  that  the  properties  of  the 
monomers  have  to  be  calculated  only  once,  after  which  the  interactions  may  be  evalu¬ 
ated  easily  and  efficiently  at  as  many  dimer  geometries  as  required.  For  this  reason  it  is 
possible  to  carry  out  much  more  accurate  calculations  than  are  possible  if  a  complete 
SCF  calculation  on  the  dimer  has  to  be  done  at  each  dimer  geometry. 


9 


MODELING  OF  ORGANOMETALLIC  REACTIVITY  USING  A  COMBINATION 
OF  EXTENDED  HUCKEL  AND  MOLECULAR  GRAPHICS  TECHNIQUES 


Pierre- Yves  Morgantini,  Peter  Fluekiger  and  Jacques  Weber 
Laboratory  of  Computational  Chemistry,  University  of  Canava, 
30  quai  Ernest  Ansermet,  1211  Geneva  4,  Switzerland 


A  new  formalism  has  been  developed  in  order  to  evaluate  the  intermolecular  interaction 
energy  between  an  organometallic  or  inorganic  substrate  S  and  an  incoming  reactant  R  in 
the  framework  of  the  extended  Huckel  (EH)  method.  Approximate  procedures  are  used  to 
estimate  electrostatic  (Ees),  charge  transfer  (E^)  and  exchange  repulsion  (Esx) 
components,  which  leads  to  short  response  times  allowing  to  use  the  model  on  an 
interactive  graphics  facility.  In  addition  to  the  electrostatic  potential  of  S.  evaluated  using 
EH  wavefunctions  and  the  NDDO  or  the  Mulliken  approximations  for  the  one-electron 
integrals,  we  are  taking  account  of  the  S-R  orbital  interactions  which  allow  to  estimate  the 
E^  component.  Finally,  the  the  short-range  Eex  exchange  repulsion  energy  is 
approximated  using  the  hard  spheres  model  for  S-R  interaction  on  the  molecular  surface 
of  S.  The  total  S-R  interaction  energy  is  then  used  as  a  reactivity  index  and  evaluated  at 
selected  points  belonging  to  the  molecular  envelope  of  S,  a  proton  with  an  empty  Is 
orbital  being  chosen  as  the  model  electrophile,  and  an  H'  hydride  ion  with  two  is 
electrons  as  the  model  nucleophile. 


Color-coded  three-dimensional  Connolly  dot  surfaces  are  used  for  the  graphics 
representation  of  the  reactivity  index  of  the  substrates,  together  with  special  procedures 
we  have  recently  developed  on  the  PS-390  in  order  to  generate  solid  models  clipped  so  as 
to  allow  the  simultaneous  visualization  of  the  structural  skeleton.  On  the  basis  cf  the 
results  obtained  for  a  large  series  of  organometallic  reactions,  this  model  is  shown  to 
describe  adequately  the  initial  stage  of  electrophilic  and  nucleophilic  addition  or 
substitution  mechanisms.  In  particular,  ferrocene  and  iron  pentacarbonyl  are  correctly 
predicted  to  undergo  an  electrophilic  attack  on  metal,  whereas  for  arene-M(CO)3  species 
and  their  derivatives,  the  nucleophilic  attack  takes  place  as  expected  on  the  exo-face  of 
the  lioand  ring.  When  applied  to  the  modellizaron  of  organometaliic  reaction  mechanisms, 
the  combination  of  simple  quantum  chemistry  methods  and  molecular  graphics 
techniques  seems  therefore  able  to  bring  an  interesting  contribution  towards  a  better 
understanding  of  the  processes  of  specific  interactions  between  chemical  species. 


GENMOL 


A  FAST  PROGRAM  FOR  MOLECULAR  MODEL I 2AT I ON 

APPLICATION  TO  THE  DETERMINATION  OF  THE  ORIGIN 
OF  THE  PSYCHOTONIC  OR  SEDATIVE  ACTIVITY  OF  TRICYCLIC 
ANTIDEPRESSANT  DRUGS  . 

Gerard  PEPS 

■1C '  -CNA3  Campus  Luminy,  case  9i3,  13233  Marseille  Cedex  3. 


Conventional  methods  of  molecular  mechanics  allow  to 
achieve  the  molecule  geometry  ny  minimization  of  tne  total 
strain  energy,  using  partial  derivatives  of  this  energy  on 
each  atom. 

-  Such  calculations  are  time  consuming  (  the  number  of 
optimization  cycles  increases  with  the  number  of  atoms  ). 

-  The  molecular  geometry  which  is  obtained  to  the  closest 
energy  minimum  of  the  starting  point,  without  any 
possibility  to  know  whether  this  minimum  is  the  deepest  one 
(  generally  che  conventional  methods  do  net  perform 
conformational  analysis  ). 

The  aim  of  our  program  named  GENMOL  is  above  all  to 
generate  the  most  stable  conformations  of  a  molecule.  The 
calculations  are  based  on  the  concept  of  pivots  which  may  be 
bonds  or  not.  A  hierarchy  of  the  deformations  is  defined  ( 
torsion,  bending  and  stretching  ).  The  program  links 
optimized  fragments  of  molecules  coming  from  differents 
libraries.  Conformational  analysis  is  performed  during  the 
building  of  the  molecule  and  the  privileged  conformations 
are  kept.  GENMOL  is  automated  for  polypeptides  and  nucleic 
acids  (  RNA  or  dna  ) . 

GENPIOL  is  one  hundred  times  faster  than  equivalent 
programs.  Bending  and  stretching  parameters  are  analoguous 
to  MM2  (  ALLINGER  )  parameters.  Torsional,  van  der  WAALS, 
hydrogen  bond  parameters  are  derived  from  ECEPP  (  SCHERAGA 
).  The  electrostatic  energy  is  computed  with  net  atomic 
charges  issued  from  DEL  RE's  method  for  the  sigma  part,  and 
with  an  empirical  approximation  for  the  PI  part  in  order  to 
get  PAP. I SER  and  PARR  values. 

When  the  molecule  conformations  of  19  antidepressant 
tricyclic  drugs  are  compared  to  those  of  the 
neurotransmitter  present  in  the  corresponding  synapses,  it 
turns  out  that  two  structurally  different  families  exist  for 
serotonin  receptors.  The  psychotonic  of  sedative  effect  are 
related  to  the  receptor  family  touched  by  the  drug. 

The  calculations  indicate  that  the  serotonin  receptors 
touched  by  the  sedative  drugs  are  structurally  similar  to 
histamine  recep' ors,  what  means  that  they  also  affect  the 
histaminic  central  nervous  system  and  contribute  to  the 
sedative  effect. 
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RN  INTRODUCTION  TO  NUMERICRL  OPTIMIZATION 

Claude  LEliRRECHRL 
INRIR  -  BP  105  -  78153  Le  Chesnay 


To  solve  an  optimization  problem  is  to  find  a  set  of  variables, 
satisfying  some  given  constraints,  and  minimizing  a  given  cost  -  function.  We 
consider  only  the  so-called  nonlinear  pr npr annina  problem,  in  which  the 
variables  vary  in  Rn,  and  there  are  finitely  many  constraints. 

The  solution  (if  there  is  one)  is  actually  approximated,  using  an 
algorithm  which  constructs  an  iterative  sequence  of  trial  solutions.  The  user 
must  provide  a  program  which,  for  any  given  value  of  the  variables,  computes 
the  corresponding  values  of  the  cost  and/or  constraints  (and  also  their 
partial  derivatives).  From  this  point  of  view,  optimization  is  fairly  similar  to 
solution  of  systems  of  nonlinear  equations. 

We  explain  the  main  ideas  that  are  used  to  define  optimization 
algorithms.  They  are  best  viewed  in  the  unconstrained  case,  i.e.  when  the 
variables  are  allowed  to  take  any  value  in  Rn.  In  this  framework,  the  simplest 
problem  is  when  the  cast-  fv  tctian  is  quadratic,  say 

f ( x )  =  2*tFIx  +  bT x  +  c  . 

Here,  fl  is  the  matrix  of  second  derivatives  of  f,  which  actually  do  not  depend 
on  x.  The  problem  then  reduces  to  a  linear  system  of  equations:  Ax  +  b  =  0. 

When  f  is  not  quadratic,  its  matrix  of  second  derivatives  varies  with  x. 
To  be  efficient,  an  optimization  algorithm  must  somehow  estimate  this  matrix; 
this  gives  a  quadratic  function,  which  can  be  used  to  approximate  f.  This  is  the 
basis  for  Newton-like  methods,  conjugate  gradient,  etc...  We  demonstrate  this 
approach  and  we  explain  how  it  can  be  adapted  in  various  situations,  such  as 
when  the  problem  is  large-scale,  or  stiff,  or  when  it  originates  from  a  least- 
squares  problem,  etc... 

Finally,  we  show  how  the  above  ideas  can  be  extended  to  constrained 
situations:  first  in  the  case  of  linear  constraints  (the  feasible  domain  is  a 
polyhedron),  and  also  in  the  general  case  when  the  constraints  are  nonlinear. 

In  the  former  case,  th?  sequence  of  iterates  is  usually  feasible,  while  it  can 
hardly  be  so  in  the  latter:  it  is  only  asymptotically  that  nonlinear  constraints 
are  eventually  satisfied. 
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Conformational  analysis  of  macrocyclic  receptors 
and  of  their  substrate  complexes. 

From  static  to  dynamic  models. 


G.  WIPFF, 

UA  422  du  CNRS 

Instirut  de  Chimie 

1,  rue  Blaise  Pascal 

67000  STRASBOURG  (France) 


Macrocyclic  synthetic  receptors  and  their  complexes  allow  to  test  cur¬ 
rent  theoretical  and  computer  graphics  techniques,  and  to  provide  a  deeper  in¬ 
sight  into  features  which  are  essential  for  molecular  recognition,  but  which 
may  not  be  easily  amenable  to  experimental  studies  :  conformational  preferen¬ 
ces  of  the  free  and  complexed  receptor,  effects  of  solvent  and  environment  on 
structure,  binding  affinity  and  selectivity,  macro(poly)cyclic  effect,  nature  and 
orientation  of  the  binding  sites,  size  of  the  concavity,  etc.... 

We  combine  molecular  mechanics,  molecular  dynamics  and  Monte 
Carlo  simulations,  normal  modes  analysis,  to  adress  the  above  questions  on 
prototype  macrocyclic  compounds  :  crown  ethers,  bicyclic  and  tricyclic  cryp- 
tands,  cryptates  and  derivatives.  From  the  collection  of  optimized  structures,  as 
well  as  from  the  dynamics  studies,  we  get  dynamic  views  of  these  systems. 
Their  flexibility  /  rigidity  is  discussed  as  a  function  of  conformation,  connecti¬ 
vity,  and  topology  of  the  free  /  complexed  state,  and  the  "Receptor-Substrate" 
complementarity. 

From  a  technical  point  of  view,  several  force  field  representations  are 
compared,  and  the  crucial  role  of  a  high  quality  computer  graphics  system  to 
build  3D  structures  of  "supermolecules",  and  to  analyze  static  and  dynamic  re¬ 
sults  is  emphasized. 


INTERMOLECl  LAR  FORCES  [>  SOME  IONIC  ATMOSPHERIC  CLUSTERS  INVOLVING  WATER 
QL  AMTM  MECHANICAL  and  MONTE  CARLO  CALCULATIONS 

E.  Kochanski  A.  Rahmouni  R.  Wiest  2 

ER  ni39  du  CNRS.  Laboratoires  do  Chimie  Theorique  •  et  Quantique.  :  . 

Universite  Louis  Pasteur.  Institut  de  Chimie.  BP  296  RS,  6-’008  Strasbourg.  France 

Quantum  mechanical  and  Monte  Carlo  calculations  huso  beer,  perfurncq  for  the 
atmospheric  ionic  clusters  HiO'{HL0)n  and  Cir\n:0)„  .  n  =  1  to  9  or  !0.  These  systems 
being  stabilized  by  interm.olecular  interactions,  special  attention  has  been  paid  to  a  correct 
description  of  the  different  energetic  contributions. 

In  a  first  step,  accurate  quantum  mechanical  calculations  (SCF  and  dispersion 
energies)  hate  been  performed  for  many  geometrical  configurations  of  bimolecuiar  systems 
M' ...H  O  I M~  =  H-O-  and  Grl  .  allowing  the  determination  of  parameters  in  order  to 
represent  the  quantum  mechanical  results  by  analytical  expressions.  Such  expressions  are 
used  in  Monte  Carlo  calculations,  to  compute  the  energy  of  any  geometry  of  the  clusters. 

In  most  of  the  Monte  Carlo  calculations  presented  in  this  work,  the  energy  of  the 
dusters  is  represented  as  a  sum  of  the  interaction  energies  between  pairs  of  molecules.  An 
'atom-atom'  (I2-I-M  formula  .  with  approriate  parameters,  has  been  used  to  describe  the 
interaction  between  the  ion  and  one  water  molecule.  Different  fittings  of  the  ab  initio 
values  arc  proposed,  allowing  to  study  the  parameters  effect  on  Monte  Carlo  results. 

The  interaction  between  two  water  molecules  is  taken  from  the  work  of 
Matsuoka  et  all  (ref.  1). 

Three-body  contributions  are  also  considered,  this  effect  being  very  important  in 
some  cases,  in  particular  to  describe  the  clustering  energies  of  very  small  clusters. 

Two  geometries  of  H}0~  are  considered  :  a  planar  and  a  slightly  pyramidal  one. 

In  the  second  step  of  the  work.  Monte  Carlo  calculations  have  been  performed 
with  the  different  sets  of  parameters  determined  in  step  1.  Energetic  and  structural  in¬ 
formation  is  then  obtained  and  analyzed.  A  special  attention  is  paid  to  clustering  energies, 
which  can  be  compared  to  experimental  values,  and  to  the  distribution  of  the  water 
molecules  around  the  ion.  The  statistical  distribution  of  the  water  molecules  and  the 
geometrical  structure  of  the  most  stable  configurations  are  considered  . 

In  a  third  step,  we  are  checking  the  validity  of  the  approximations  by  performing 
ab  initio  calculations  on  some  geometrical  structures  of  special  interest  generated  in  the 
course  of  Monte  Carlo  calculations.  SCF  calculations  are  now  performed  on  the  whole 
clusters  and  not  only  on  the  bi  or  trimolecular  systems.  The  nature  of  the  filling  of  the  first 
soivation  shell  deserves  a  particular  attention.  Let  us  note  that  we  are  interested  in 
particular  cases  and  not  only  in  average  values.  Comparison  with  other  parameters  or  ex¬ 
pression:,  proposed  in  the  literature  will  also  be  done. 

Part  of  this  work  has  already  been  published  recently  (ref.  2). 

(1)0.  Matsuoka.  M.  Yoshimine,  E.  Clementi,  J.  Chem.  Phvs.,1976,  64  ,1351 
S.  Romano.  E.  Clementi,  Gazz.  Chim.  Ital.,1978.  J_08  ,319 
(  2)  E.  Kochanski,  E.  Constantin,  J.  Chem.  Phys.,  1987,  87  .1661; 

E.  Kochanski,  Nouv.  J.  Chimie,  1984.  8  .605;  Chem.  Phvs.  Letters,  1986.  L30  .291: 
J.  Amer.  Chem.  Soc.,  1985,  J07  ,7869;  Chem.  Phys.  Letters.  1987,  JT3  ,143 
E.  Kochanski,  A.  Rahmouni.  J.  Chim.  Phys.,  on  press 
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Etude  comparative  des  interactions  intermolSculaires 
au  noyen  de  m€thodes  pertubatives  et  variationnelles. 


Application  au  dim&re  de  l’eau. 


O.Hess,  J.Caillet,  J.Langlet 

Dynamique  des  Interactions  Moiacuiaires 
Universita  Pierre  et  Marie  Curie 


Ce  travail  concerne  la  theorie  des  interactions 
mtermoieculaires  au  niveau  iondamentai ,  utiiisant  des 
methodes  ab  initio  elaborees,  aussi  bien  de  type  perturbatif 
que  variationnei  avec  utilisation  de  bases  d’orbitales  ato- 
miques  de  taille  tres  importante. 

Au  niveau  perturbatif , nous  avons  mis  en  evidence  le  role 
important  des  contributions  d'echange  du  second  ordre  ( d ’ in¬ 
duction  et  de  dispersion)  qui  sont  habituellement  negligees 
dans  les  calculs  d ' interacti ons  intermoleculaires .  En  parti- 
culier,  il  faut  done  egalement  tenir  compte  de  la  dispersion 
d'echange  lorsque  1 '  on  rajoute  les  contributions  de 
correlation  intermcleculaire  a  1'energie  d' interaction. 

Cependant  la  comparaison  avec  les  calculs  de  type 
variationnei  montre  que  dans  une  approche  purement  pertur¬ 
bative  les  contributions  d' ordre  superieur  jouent  un  role  non 
negligeable,  meme  dans  les  regions  d'equilibre  du  complexe. 

P  ^  zlz.  2. 1  c. "11“ c ,  ”  r 9 s ? r. t c ^ r  o  L  cp  i  ^  ^  r  9  9'^  1 1 9 1 =  vi  v*  ~ 
trant  la  dependance  des  differentes  composantes  de  1'energie 
d ' interaction  avec  la  geometrie  du  dimere  (differentes 
orientations  relatives). 

De  nouveaux  potentiels  semi-empiriques  calibres  a  partir 
de  cette  etude  sont  proposes. 
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MODEL  HAMILTONIAN  AND  PERIODIC  HARTREE-FOCK 
CALCULATIONS: 

WATER  molecule  in  different  surroundings 


J.  G.  ANGYAN  and  3.  SILVI 

Laboratoire  de  Spectrockimie  \foUcuiaire,  f ERA  508),  Vniversite  Pierre  et  Mane 
Curie,  Place  Jussieu.  7 5230  PARIS  Cedes,  ( France ) 


ABSTRACT 

Several  approaches  can  be  used  to  describe  the  water  molecule  in  dif¬ 
ferent  physical  and  chemical  environments.  On  the  one  hand  is  the  peri¬ 
odic  Hartree-Fock  method  [l;  which  applies  to  crystal  phase  cases  such  as 
proton-ordered  ice  and  LiOH.HoO  providing  ab  initio  reference  wavefunc- 
tions.  On  the  other  hand  are  model  (or  effective)  hamiltonian  methods  in 
which  the  surrounding  is  represented  by  an  effective  potential.  The  self- 
consistent  Madelung  potential  (SCMP)  method  ,2]  uses  a  specific  form  of 
the  effective  potential,  involving  lattice  sums,  adapted  to  ordered  systems. 
Introduction  of  a  general  form  of  the  atom-atom  pair  distribution  function 
allows  to  extend  this  latter  model  to  disordered  situations  such  as  pure 
molecular  liquids.  A  series  of  calculations  of  the  water  molecule  in  different 
environments  has  been  performed.  Charge  distributions  (atomic  multipoles 
and  electron  density  maps)  computed  by  the  two  methods  are  compared  and 
analyzed. 

REFERENCES 

1]  C.  PISANI,  R.  DOVESI  and  C.  ROETTI,  Lecture  notes  in  Chemistry,  48, 
(1988) 

'2j  J.G.  ANGYAN  and  B.  SILVI,  J.  Chem  Phys.,  86,  6957  (1987) 


STATISTICAL  MECHANICS  APPROACH  IN  MCBSLjNG  -OlECULAR  jTRUCTUR  £S 
AND  CONFORMATION -DEPENDENT  PROPERTIES 

S.  S cordama g l i a ,  L.  Bar  ino 

Himont  I  tat  ia  i  Centro  S  ictrche  Ncvara,  Via  Caduti  del  Lavoro, 

23100  NOVARA  (Italy) 

ABSTRACT 

An  averaged  picture  of  the  conformational  possibilities  of  a 
molecule  while  approaching  its  own  interactive  site  gan  be 
calculated  through  the  statistical  analysis  of  all  the  acceptaole 
molecular  conformations . 

The  CSD  (Conformations  Statistical  Distribution)  method  is 
presented  which  calculates  the  whole  rotational  hypersurface  of  a 
molecule  having  n  internal  rotational  degrees  of  freedom,  through 
the  calculation  of  the  nonboncea  intranolecular  energy 

( Lenn a r d- J on e s  term  plus  coulombic  interaction  between  nonbonded 
atoms)  for  each  point  of  a  n *d i mens i o na l  homogeneous  grid  in  the 
surface,  and  then  applies  Boltzmann's  statistical  mechanics  to  all 
the  allowed  co nf ormat i ona l  microstates  thus  obtaining  the 
statistical  weights  of  the  various  minima  in  which  the  molecule  can 
be  found,  the  partition  function,  entropy  and  thermodynamic 
potentials  both  in  the  hypersurface  and  around  the  detected 
n ini ma . 

Some  examples  shew  the  use  of  the  calculated  distribution  functions 
in  the  determination  of  some  weighted  con f o rm a t i on -dep enae n t 
molecular  properties  such  as  internal  rotational  freedom,  molecular 
shape,  total  dipole  moment  and  other  steric  ana  electronic 
features.  These  cuantities  appear  suited  molecular  descriptors  in 
many  theoretical  studies  such  as  prediction  of  most  probable  3-D 
structure  in  isolated  mo  l  e c u  le s , s t r u c t ur e-a c t i v i t y  ana 

s t rue t ure-func t i cn  relationships,  modeling  of  molecular 
interactions  with  a  solid  surface,  evaluation  of  intramolecular 
physical  properties  affecting  packing  between  molecules  and  its 
modifications  with  temperature.  In  particular  for  representative 
models  of  flexible  chain  molecules  the  effects  linked  to  the 
t rans f orma t i on  of  thermal  energy  absorbed  by  a  compound  into 
internal  rotational  energy  can  be  analyzed  pointing  out  chain 


c onf o rma t i ona l  transitions  versus  t e« pe ra t u r e  .  The se  are  correlated 
to  phase  transitions  in  the  c o r r e s pond ing  polymeric  materials. 


ASPECTS  DB  LA  DYNAMIQUB  DBS  IONS  EN  SOLUTION  PAR 
SPBCTROSCOPIB  IR  LOINTAIN 

B  GUILLOT  ,  Ph.  MARTEAU  et  J.  OBRIOT 

Laboratoire  de  Physique  Theorique  des  Liquides. 

Universite  Pierre  et  Marie  Curie. 

4  Place  Jussieu,  Paris  75005. 

Laboratoire  des  Interactions  Moleculaires  et  des  Hautes  Pressions. 

CNRS  Universite  Paris-Nord,  Avenue  J.B  Clement. 

Villetaneuse  93430 

Les  spectres  infrarouge  lointain  des  halogenures  alcalins 
Lid  ,  Lil,  Nal ,  Csl  et  Rbl  en  solution  dans  1 '  acetoni  t  ri  le ,  l'acetone  et 
le  methanol  ont  ete  enregistres  dans  le  domaine  spectral  30-600cm  .  Ces 

spectres  presentent  une  forme  de  bande  complexe  caractensee,  en 
general,  par  plusieurs  maxima  d'absorption  (ex.  pour  HI/'CH30H  trois 
maxima  a  35,  235,  460  cm  1  ) 

Une  analyse  theorique  permet  d 'exprimer  la  densite 
spectrale  a  l'aide  (i)  de  la  fonction  de  correlation  du  champ 

polarisant,  emanant  de  la  charge  ionique,  lorsque  les  ions  sont 
dissocies,  (ii)  de  la  fonction  de  correlation  des  moments  dipolaires  des 
agregats  ioniques  non  dissocies  et  (iii)  d ' une  fonction  de  relaxation 
structurale  caracteri sant  le  rearrangement  des  molecules  de  solvant 
autour  de  1 1  ion  La  fonction  de  correlation  du  champ  polarisant  est 
ensuite  reduite  a  une  expression  plus  manipulable  en  postulant  que  le 
mouveraent  relatif  des  ions  par  rapport  aux  molecules  de  solvant  est 
caracterise,  aux  temps  courts,  par  une  oscillation  dans  la  cage  des 
premiers  voisins  et  aux  temps  longs.  par  un  mouvement  de  diffusion 
t  rans la t  lomiel  L 1  appl  ica  t  ion  de  la  theorie  de  Zwanc ig-Mor i  a  la 

fonction  de  correlation  du  champ  polarisant  ainsi  qua  celle  des  ions 
non  dissocies  permet  de  rendre  compte  quantitat ivement  des  spectres  IRL 
des  ions  en  solution.  En  particulier,  en  utilisant  le  concept  de 
potentiel  de  force  raoyenne,  il  est  possible  d'attribuer  aux  maxima 
d'absorption  les  frequences  de  vibration  des  anions  et  des  cations  dans 
leur  cage  respective,  d'une  part,  et  celles  des  paires  d'lons  en  contact 
ou  separes  par  une  molecule  de  solvant,  d 'autre  part.  Ce  dernier 
mecanisme  est  decrit  comme  le  resultat  dun  echange  lent  entre  deux 

etats  energetiques  separes  par  une  barriers  de  potentiel  de  quelques  kf 


Determination  of  nucleic  acid  conformations 
by  vibrational  spectroscopy  accompanied  with 

harmonic  dynamics  calculations 


M.  Ghoai ,  R.  Letellier,  D.  Dohy  &  E.  Taillandier 

Laboratoire  da  Spectroscopic  Biomoleculaire. 

U.F.R.  Biomedicale  da  Bobigny.  Universite  Paris  XIII 
74,  rua  Marcal  Cachin.  93012  Bobigny  cedes.  Franca 


It  is  now  wall  known  that  tha  vibrational  spectroscopy  can  be 
considered  as  a  powerful  probe  to  identify  tha  secondary  structure  or 
oligo-  and  polynucleotides.  In  addition  of  tha  experimental  spectra 
interpretation,  the  normal  coordinate  analysis  allows  us  to  dive  a 
quantitative  evaluation  of  nucleic  acid  structural  parameters. 

In  our  previous  studies,  we  have  carried  out  a  systematic 
investigation  on  DNA  vibration  modes  (1).  The  Wilson's  GF  method  has  been 
applied  on  the  dynamic  models  constituted  by  nucleosidic  residues.  In 
this  approach,  a  non- redundant  valence  force  field  has  been  employed. 
Recently,  these  calculations  have  been  extended  in  order  to  study  the  V- 
deoxyribose  modes  as  a  function  of  the  pucker  pseudorotational  parameters 
(2)  . 

By  the  use  of  this  reliable  force  field  through  the  Higgs  method 
(infinite  polymeric  helical  chain  phonon  dispersion  curves),  we  have  also 
performed  a  critical  discussion  on  the  conformations  adopted  by  poly 
idG) ,  poly  (dC)  in  solution  at  various  salt  concentration  (3). 

The  present  work  is  concerned  with  the  quantitative  evaluation 
of  the  sugar  conformation  (P  angle) ,  as  well  as  the  purine  base 
orientation  (  “X  angle)  .  The  discussion  will  be  carried  out  on  the  manner 
how  the  (P,%)  couple  of  structural  parameters  can  be  determined  in  both 
eC  and  B  deoxynucleosides  by  the  use  of  the  vibrational  markers. 
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FORCE  FIELD  CALCULATIONS  OF  SIDEROPHORES. 

STABILITY  AND  CONFORMATIONS  OF  Fe111  CHELATES  WITH  NOVEL  IRON 
RELATED  PARAMETERS. 


A.3ouraouia.  M.Fcthailan3.  F.Mhennm,.B.3lcivea  ana  R.Gaiio0. 
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Iron  is  essential  in  human  where  an  excess  is  toxic,  and  in  plants  where  a  lack  of 
the  element  causes  ferric  chlorosis. 

Iron  carriers  (siderophores)  are  appropriate  organic  ligands  aple  to  complex  and 
transport  Fe111  for  establishing  an  appropriate  concentration  level  of  the  element.  We  have 
therefore  undertaken  the  design  and  synthesis  of  new  efficient  iron  cheiatants. 

The  new  ligands  contain  carboxy  and  catechol  substructures  arranged  in  order 
that  six  oxygen  atoms  are  located  at  the  top  of  an  octahedral  structure.  The  use  of  a 
"molecular  model"  approach,  for  the  choice  of  the  organic  framework  of  the  ligand, 
ceing  an  oversimplification,  we  have  made  a  detailed  force  field  calculation  analysis  with 
the  Molecular  Mechanics  program  BLEMO  using  the  1985  parametrization  of  the  ALLINGER 
MM2.  These  calculations  have  needed  the  design  and  testing  of  new  parameters,  net  yet 
available,  for  describing  the  interactions  with  iron. 

The  study  of  the  iron  complexes  as  a  function  of  the  ligand  size  and  of  the  nature 


of  the  possible  isomers  shows  what  are  the  most  stable  structures,  and  helps  in  tailoring 
other  ligands.  The  syntheses  of  these  new  siderophores  are  in  the  course. 


Handling  the  structural  information  arising  from 
molecular  modeling  for  drug  activity  predictions  : 
the  SARAH  system. 
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The  computer  is  a  versatile  tool.  In  particular  it  can  perform  either  numerical  compu¬ 
tations  or  logical  operations  leading  to  a  branch  of  Artificial  Intelligence.  In  the  case  of 
molecules  like  drugs  acting  in  complex  processes,  both  of  these  capabilities  may  be 
used  to  set  up  a  computer  assisted  drug  design  system.  The  molecular  properties  : 
geometry,  conformational  energies  and  also  electronic  properties  can  be  computed  by 
means  of  molecular  mechanics  and  quantum  chemistry,  but  handling  this  very  rich 
information  efficiently  requires  special  procedures  which  may  mimick  the  chemist's 
reasoning  with  the  advantages  of  getting  rid  of  subjectivity  and  performing  systematic 
and  lenghthy  comparisons  without  the  risks  of  human  fatigue. 

The  SARAH  system  (Structure-Activity  Relationships  by  Apprenticeship  and 
Heuristics)  rests  upon  the  Key-Lock  model  :  the  substrate  has  to  fit  the  receptor's  site, 
in  particular  through  a  special  molecular  pattern  and  a  shape  adapted  to  the 
receptor's  geometry.  It  also  has  to  exhibit  the  proper  electronic  properties  to  make  the 
substrate-receptor  complex  as  active  as  possible. 

The  procedure  is  decomposed  into  three  steps  : 

1.  A  characteristic  pattern  of  heteroatoms  is  defined  by  superimposition  of  the  active 
molecular  skeletons  to  evidence  the  recurrent  atoms. 

2.  Rules  of  activity  are  determined  by  systematically  checking  the  correlations 
between  the  variations  of  shape  or  electronic  properties  around  the  characteristic 
pattern  and  the  variations  of  activity  relative  to  a  reference  molecule. 

3.  A  multivariate  statistical  analysis  is  performed  on  the  electronic  properties  of 
active  molecules  to  extract  a  quantitative  criterion  for  the  activity  of  a  new 
molecule. 

The  system  is  evolutive  :  when  a  new  molecule  has  been  tested,  it  can  enter  the 
training  set  and  refined  criteria  are  produced. 


The  possibilities  of  the  system  are  illustrated  by  the  analysis  of  the  antiepileptic 
properties  of  a  series  of  benzodiazepines. 
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RELATION  ENTRE  CONSTANTES  DE  FORCE  ET  CONSTANTES 
ELASTIQUES  DANS  DIVERS  MATERIAUX 
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Le  but  de  ce  travail  est  de  relier  1’elasticiie  mieroseopique  (celle  des  liaisons)  et 
i'elasticite  macroscopique.  Une  methode  de  caicui  tnatriciel  des  elements  Cn  du  tenseur 
d’elasticite  a  partir  des  constantes  de  force  a  etc  developpee  par  Shimanouchi  en  1971. 
reposant  sur  la  relation  generate  de  Bom  et  Huang. 

Nous  avons  ecrit  un  programme  de  calcul,  a  partir  de  cette  theorie,  permettant  de 
calculer  les  elements  Cjj  4  partir  du  champ  de  force  d’un  cristal.  Ce  calcul  donne 
egalement  les  derivees  partielles  des  constantes  elastiques  par  rapport  aux  constantes  de 
force,  c’est-a-dire  la  contribution  de  chaque  type  de  liaison  et  de  chaque  type  d’angle 
dans  la  valeur  de  la  constante  elastique. 

Nous  avons  etudie  des  materiaux  trfes  divers  dont  les  constantes  elastiques  sonr 
connues  et  les  resultats  ont  montre  que  : 

1°)  les  champs  de  valence  utilises  sont  tout  4  fait  vaiables  pour  calculer  avec  une  bonne 
approximation  les  constantes  elastiques  et  les  frequences  de  vibration. 

2C)  la  contribution  relative  de  chaque  liaison  et  angle  permet  de  mettre  en  evidence  les 
param£tres  qui  influencent  le  plus  le  comportement  elastique  des  materiaux  : 

a)  les  resultats  obtenus  pour  le  diamant,  le  graphite  et  le  polyethylene,  montrent  que 
les  Cjj  sont  tr4s  sensibles  4  la  dimensionality  du  reseau  covalent. 

b)  les  resultats  compares  du  quartz  alpha  et  de  l’alumine  alpha  montrent  I'influence 
de  la  coordinence  des  atomes. 

c)  les  resultats  obtenus  pour  BaTiO^  montrent  I'influence  de  l'orientation  des 
liaisons  par  rapport  aux  axes  du  cristal.  Lorsque  les  liaisons  sont  paralleles  aux  axes  du 
cristal,  les  constantes  elastiques  dependent  pratiquement  exclusivement  de  celles-ci,  du 
moins  en  ce  qui  concerne  la  compressibility. 
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5  IMITATION  OF  FREE  ENERGY  RE  L  AT  I  ONSH  I  P  5  FOR  THE  ELECTRON  TR\N-iFH< 
REACTION  BETWEEN  Fe~  and  Fe2~. 

A.  Gon/al  c/  -  Lar’onr  .  I.M.  1  Inch.  \.  OLiva  and  /.  3«m  r  :  m  . 

Depu  it  rm?n  t  do  (juunica.  Ini  v-t  -»  i  *  at  \nr onrima  d<*  3.t  i  u  »_*  I  ■  •  n.i  .  ■'"I 

Be  L  La  ro  r  r  a  i  Barcelona)  ,  Spain 

The  classical  treatment,  or'  electron  tran-sTer  :  <‘ai/  i  on 
•lenerallv  use  an  act  i  vateci-eomp  Lex  rarntalisra.  \ccordi  nj:  r  o  r  h«* 
F  ranck -Condon  principle,  the  electron  transfer  musr  >«<;ur  i*  »!'.*■ 
intersection  region  of  the  diabaUc  putcin.  ial  c.i v,:  u;. 
,‘iyper -surfaces .  This  region  is  reached  by  a  suitable  fluctuar  ion 
m  the  nuclear  configurations  of  the  precursor  complex. 

The  expression  of  the  rare  constant  given  by  the 

standard  r  ran  s  i  t  i  on  -  t  a  t  e  rheorv  shows  that  it  depends  on  the 

exponent,  ial  fact  or  exp  i  -  ^  •*  >  which  mejiun*.  t.h<»  probab  i  li  c  v 

that  the  system  will  bo  in  the  r  rans i t  Lon-state  region  relative 
r  .i  the  probability  of  being  at  the  reactants,  region.  The 
ova  lu.it  Lon  of  the  react  ion  rare  would  requite  an  extremely  Lons 
iimul.ti  Lon  time  if  one  started  at  r  lie  react  jm  region. 

In  this  work  we  present  the  results  obtained  for  fit’ 
electron  transfer  reaction  between  Fe  and  Fe“  .  rite  distance 
between  both  ions  bema  kept:  fixed.  Given  that  the  two  reactants 
ite  mo  non  t  oin  l  c  ions.  the  nuclear  coordinates  of  t  h  ;  sxstem  on  1  s 
correspond  to  solvent  coordinates.  For  a  sample  containing  i 0 
witer  mo  )  .  -  e  u  1  •  •  s  we  h  a  \  a.  -  no :  i 1  e  (1  one  million  oi  enfum.c  min 

i '  a :  t  i  n  p  it  .  >m  the  tear’. nits  t>'iMn.  hut  none  .it'  '  hem  has  be.  •  n 

round  a '  the  transition  stare  n'j  i  .in. 

To  e v  a  1 ua r  e  the  free  energy  difference  we  have  used  two 
different  techniques.  Trt  one  of  them  we  have  assumed  a  parabolic 
behavior  of  the  free  energy  surface  with  respect  to  the  solvent 
reaction  coordinate.  Th  i  s  has  permitted  to  obtain  the  probability 
of  tindins  conf  i  aur.it  i  tins  at  the  transit  ion-state  region  from  the 
calculation  of  several  values  of  the  logarithm  of  the  probability 
corresponding  to  configurations  Lying  in  the  reactants  region. 
The  second  technique  has  allowed  to  obtain  a  more  quantitative 
value  using  a  free-energy  perturbat ion  technique  with  a  suitable 
mapp t  ng  pot  erit  i  a  1  . 


POLYMORPHISM  AND  MOLECULAR  DISPLACEMENTS  IN  ORGANIC  SOLID  STATE. 
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Organic  molecules  often  display  several  crystalline  varieties.  The  relative 
stability  of  such  arrangements  ^nd  their  mutual  transformation  have  been 
the  object  of  numerous  static  and  dynamic  models.  Almost  all  studies  are 
based  on  atom-atom  potential  methods.  Lattice  and  molecular  dynamic  methods 
have  been  applied  to  plastic  phase  and  to  the  corresponding  ordered- 
disordered  transition  of  small  molecules  or  rigid  molecules,  such  as 
benzene  or  naphthalene.  Relative  stability  calculations  have  been  often 
carried  out  for  testing  the  efficiencies  of  atom-atom  potential  methods. 
Dif f erencies ,  often  of  the  order  of  few  Kcal/mole,  are  generally  well 
reproduced,  even  if  absolute  values  depend  of  the  choice  of  parameters. 

Concerning  the  transition  between  phases,  it  is  to  be  said  that  experimen¬ 
tal  studies  show  a  great  variety  of  parameter  dependancy,  feeding  contro¬ 
versies.  Single  crystal-single  crystal  second  order  transformations 
sometimes  are  claimed  to  not  exist.  Nevertheless,  such  transformations 
would  be  of  interest  from  both  theoretical  and  experimental  physic  point  of 
view,  since  they  are  conceived  as  cooperative  displacements  of  molecules. 

Studies  of  relative  stabilities  and  mutual  transformaations  of  phases  are 
of  special  interest  for  ph. rmaceutical  and  industries  of  explosives. 
Polymorphism  is  a  very  usefull  support  for  modelizing  molecular  motion  in 
space.  Modelization  of  phase  transition  of  para-dichloro  benzene  and 
diacetamide  were  investigated  in  our  laboratory  and  the  entire  transition 
path  will  serve  as  typical  examples. 
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NONARCHIMEDEAN  MODEL  IMG  FOP  3 1 CS  VS  7EMS  EVOLVING  ON  M-W  HME  SCALES 
O.Iordache^  a"d  a.T.Franacool~ 

*  -cl y-.acrm'c  Insti  dure  of  Bucharest,  Ceoartnsnt  of  Chemical  Ir.cinaer^-', 
i,  Polizu  Street,  R -73 1 25  Sucnarest  12,  Romania 

'Institute  of  Physics  and  Nuclear  Engineering,  °.0.3ox  MG-6, 

R-76900  Mag urel e-Bucharest,  Romania 

The  nonarcnimedean  (n.a.)  analysis  was  used  to  describe  relaxation  orocesses 
in  drug-membrane  interactions,  sef  1 -exci  tation  of  interfaces,  the  balance  of 
entroDy  in  biosystems  which  exhibit  more  than  one  time  scale  in  their  evolution 
[1-51.  'he  method  relies  uocn  the  reel  ace-ent  of  the  real  time  t  and  :f  the  evo¬ 
lution  variable  y(t)  by  elements  of  an  n.a.  field  [1,5, 71 

=  t  +ct+  ...  -:'’t  and  Y(T)  =  y  (t)  +  cy,(t)  +.  ..  -  :'v..,(t) 

respectively,  where  -t  is  an  exoansi  on  parameter  and  y  It!  are  the  comoonents  of 

m  ■ 

Y ( T)  corresponding  to  the  M+l  tine  scales.  The  model  eos.  are  then  reolaced  bv 
their  n.a.  corresoondents  to  the  result  of  an  enriched  mathematical  structure  of 
the  evolution  variable.  Particularly,  it  is  shown  that  the  formalism  can  easily 
accommodate  nonexponential  relaxations  in  multi -comDar  tnental  systems  (e.g.  the 
action  potential  of  frog  sciatic  nerve  in  the  presence  of  anesthetics  [21,  cross- 
linking  aldehydes  [ 3l )  and  various  regimes  of  chemically  driven  interracial  tur¬ 
bulence  in  the  presence  of  unspecific  drugs  [4],  This  new  procedure  aoolied  to 
tre  oroblem  of  time  variation  of  the  enthreoy  production  leads  to  a  hierarchy  ;f 
enthrooy  balances  and  of  stability  criteria,  which  can  be  useful  in  describino 
embr yogenes i s,  nerve  excitation  and  other  cooperative  and  far  from  ecui  1  i bri urn 
phenomena,  'he  existence  of  mere  than  one  time  scale  checked  by  the  n.a.  analy¬ 
sis  suggests  that  the  excitation  or  inhibition  of  some  scale  of  time  could  be 
eventually  used  for  controlling  the  whole  process. 
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CLASSIFICATION  DE  MOLECULES  CHIMIQOES  :  PRESENTATION  D'UNE  MAOUETTE 

D'UN  SYSTEME  EXPERT 
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lar.s  Is  cacre  : '  une  ;r. : :  :a*. .  :n  aux  ssthcass  :  a.  azzra*.  :r. 
i  ur.  susteae  expert .  nous  cresentons  un  prograaae  zeraettar.t  ;e 
Ziass:  ftar  ass  innic i teurs  as  ia  pnctosuntnese  s:eae 

1  1  ,  ’  =  ;  ;  nh :  b  :  z eurs  sent  intrcduits  par  one  cestr .  1 1  :cr. 
s : so : : f : es  de  ieur  structure  aciacuiaire  'presence  de  certains 
crcuces  chiaiouss  et  suost ; tuants .  pcs:t::rs  relatives  :e  csux~ 

.  Css  ccapisaents  as  description  peuvent,  s .  sesem.  etre 
deaar.das  par  la  susteae.  I;  an  as*,  deouit  tine  prediction  ae 
.  activite  .r.hibitr.ze  a  part;:  de  regies  czr.cernar. t  :ss 
i.-.fprsat  esae.nt :  ei  1  eaer. c  ::rc  i  eg; cues .  La  susteae  reut 

ecaieaent  tourmr  des  sense  i  cneaents  coup  iesen ta ;  res  detenus 
car  ia  case  de  connaissances  bi bi icgrapnie  par  exeapie  . 

ls  generateur  de  susteae  expert  utilise  cans  cs  t.-avai. 
est  . ’puts |  SPIRAL  develcppe  au  LEA  Matheaa: iques  Appl:quees 
oar  i .  3CUCKET .  Cast  un  .angaae  dec.aratit  inspire  de  rF.Ci.iG 
avec  le  foraaiisae  de  ia  logtque  bu  preaier  rrcre.  :Jr. 
f:.*aa!  :sae  c b jet  u  a  ate  ir.trccuit:  cn  peut  u  defir.ir  des 
classes  pretotupes;  et  des  instances  ds  classes  ■  afcjets  peur 
acdeiiser  Is  aende  etudie  qu;  est  aan -pule  car  .as  regies 
d' inference,  ia  version  pour  aicrc-ordmateurs  ccapatibles  ru  a 
ate  surtcut  eacicuee.  Lee  scceis  s  des  rcncticns  exterr.es  cr.r 
ete  utilises  pour  la  gestian  ce  1  scran  at  .  aft  ichage  ua 
tax  tes  . 
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SYSTEME  CCTANG  :  UN  DESCRUTEUR  3D  POLYVALENT 
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La  reconnaissance  en  3D  de  fragments  structuraux  reels 
ou  standard,  ou  de  femes  electroniques  (Grbi tales  Moiecuiaires 
Potential  Slectrostatique . . . )  est  &  la  base  de  r.cmbreuses 
strategies  de  simulation  molacuiaire  :  interrogation 
structurala,  simulation  spectraia,  recherche  de  pharmacophore 
et  Drug  Design. 


Nous  etudiens  actueiisnent  un  system  de  description 
des  structures  moiecuiaires  qui  repose  sur  ur.e  specification 
qualitative  des  positions  relatives  des  atones.  Ainsi  plutot 
que  de  ccmparer  les  coordcnnees  cartes iennes  de  deux  points  A 
et  B,  on  precise  que  B  est  par  exemple  situe  au  dessus  et  4 
droite  de  A.  Pour  cela,  nous  utilisens  les  concepts  d' octant  et 
de  graphe  chimique  :  on  '/alue  un  arc  i-j  non  plus  par  la 
distance  interatanique  mats  par  le  numero  de  1' octant  localise 
stir  i  et  contenant  j .  Associi  A  une  description  en  chemins  des 
graphes  chimiques ,  ce  system  dit  "0CIAN3"  foumit  un 
descripteur  de  forme  3D  qui  optimise  les  manipulations 
informatiques  et  permet  la  reconnaissance  structurale  sur  de 
•/dri tables  cri tires  de  geonetrie  tridimentionnelle . 


Le  systime  "CCTANG"  est  actuellanent  implements  dans  le 
cadre  du  Drug  Design.  TOutefois  cette  apprcche  novatrice, 
fondSe  sur  une  perception  de  i'envircnnemar.t  local  des  atones 
doit  offrir  des  resssources  importantes  dans  le  trai tenant  de 
tous  les  probliknes  relevant  de  1 '  identification  et  de  la 
gestion  de  1 ' inf ormaticn  structurale  tridimensionnelle. 


DESIGN  OF  NEW  EFFICIENT  RADICAL  CLOCKS  WITH  THE  HELP  OF 
MOLECULAR  MECHANICS 


J.-M.  MATTALIA.  M.  FATHALLAH.  A.  SAMAT.  B.  BLAIVE 
Faculte  lies  Sciences  de  Sc  Jerome.  U.A.  126  CSRS.  13397  Marseille  Cedex  13. 


'"Radical  clocks  1"  names  molecules  able  to  trap  by  intramolecular  cvclizadon  a  radical 
generated  in  an  adequate  position.  We  have  synthesized  these  last  years"  2-norbomenes  which  are 
very  efficient  a-sulfonyl  radical  traps,  i.e.  with  high  rates  of  intramolecular  addition  of  the  radical  on 
the  double  bond  fRCyCi=  -•  lO^s'^  and  lO^s'^at  80CC  for  1  and  2  respectively).  X-ray  analysis  nas 
shown  that  the  S-fe.toj-phenyl  substituent  decreases  the  C4C5S  angle  and  consequently  the  minimal 
distance  available  between  the  reactive  sites  (2.50A  (CgH^vs  2.70A(CN)).  It  seems  that  the  rate  of 
cyclization  is  correlated  to  this  minimal  distance. 


so?c;ch3)2 

1 :  R  -  CN 
2  :  R  -  C6H5 


X  =  S02,  C-O,  CH2 
Rs=H.  B 

R7 .  R'7  -  H,  Me.  iPr,  -(CH2)„-  n:2,4.5 


In  order  to  design  new  radical  clocks  with  greatest  efficiency,  we  performed  molecular 
mechanic  calculations  on  a  series  of  molecules  related  to  the  norbomene  one  (  A  ,  B  and  few  others) 
and  studied  the  influence  of  the  substitution  on  the  distance  between  the  reactive  sites.  We  used  the 
program  BLEMO  in  the  version  including  the  MM2  parametrisation-v  Calculations  were  not 
performed  on  the  radicals  but  on  the  parent  compounds  and  the  minimum  distance  (dj.q)  was 
evaluated,  after  minimization  ot  the  energy,  by  locking  the  C},  Cj,  Xg  and  Cg  atoms  in  the  same 
plane. 


Data  show  the  importance  of  the  values  of  the  C3C4C5  and  C4CjXg  angles  on  d^.g  :  a 
weak  tightening  leads  to  an  important  decrease  of  the  distance.  For  this  purpose,  substitution  on  C5 
and  C7  ( A )  and  dimethanonaphtalene  (  ]£  )  are  the  most  efficient.  This  kind  of  study  involving  a 
high  precision  about  the  calculated  structures  is  perfectly  ensured  by  the  molecular  mechanics. 
References: 
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Reconnaisance  tridimensionnelle  de  pharmacophores. 
Application  de  la  m&hode  k  la  comparaison  morphine-enkephaiine. 


L.  Moy  et  C.  Riche. 

Laboratoire  de  Cristallographie,  I.C.S.N  au  C.N.R.S  de  Gif /Yvette,  91190. 


Par  le  biais  d'une  technique  informatique,  une  methode  de  reconnaissance  tridimensionnelle  a 
ete  developpee  et  appliquee  k  des  composes  de  meme  activite  pour  lesquels  la  topologie  du  recepteur 
n'a  pa*  ete  elucidee. 

Essentiellement  basee  sur  des  criteres  geometriques,  la  methode  urilise  une  nouvelle  definition 
de  la  molecule  et  de  ses  groupements  fonctionnels  afin  de  prendre  en  compte  les  possibles 
interactions  avec  le  recepteur.  Deux  types  de  centres  d'interaction  (centres  actifs)  permettent  de 
coder  chaque  molecule  : 

-  les  atomes  eux-memes  qui  sont  utilises  dans  la  plupart  des  comparaisons. 

-  les  centres  actifs  complementaires  localises  sur  le  recepteur  (paire  electronique  d'un 
heteroatome,  centre  et  direction  d'un  phenyle). 

En  considerant  une  molecule  rigide  et  son  analogue  flexible  actif,  le  principe  de  la  methode 
consiste  k  comparer  deux  ensembles  de  centres  actifs,  un  ensemble  etant  fixe  (molecule  rigide), 
iautre  etant  renouvele  pour  chaque  conformation  de  l'analogue  flexible.  Afin  d'eliminer  le  probleme 
de  l'origine  des  coordonnees,  le  produit  d’autoconvolution  permet  de  transformer  deux  centres 
actifs  A  et  B  d'origine  O  en  un  vecteur  AB  caractdrise  par  une  norme  et  une  direction.  D'une  fa?on 
generate,  k  un  ensemble  de  N  centres  actifs  est  associe  une  liste  de  N.(N-l)/2  vecteurs,  renouvelee 
pour  chaque  nouvelle  conformation  de  l'analogue  flexible.  Les  differents  vecteurs  sont  alors  testes 
deux  &  deux  au  niveau  de  leur  produit  scalaire  (distance  entre  centres  actifs)  et  de  leur  produit 
vectoriel  (orientation  des  centres  actifs  dans  l'espace). 

L'introduction  d'un  critere  energetique  base  sur  les  interactions  entre  atomes  non  lies  ainsi  que 
le  balayage  de  l'espace  par  la  technique  des  nombres  magiques  (derivee  des  methodes  directes  en 
cristallographie)  permettent  de  reduire  considerablement  le  nombre  de  solutions  proposees  par  le 
programme.  L'originalite  du  logiciel  PPSP3  reside  dans  le  fait  qu'il  autorise  la  superposition  d'un 
nombre  iilimite  de  centres  actifs. 

L'application  de  la  methode  a  la  comparaison  morphine  -  Leu  enkephaline  fait  clairement 
ressortir  deux  hypotheses  de  pharmacophores.  L'analvse  conformationnelle  du  pentapeptide.  basee 
sur  des  criteres  energetiques  et  les  correlations  avec  sa  structure  cristalline  (Pr.  Marot,  Nancy. 
1988)  permettent  de  comparer  les  deux  solutions.  Les  valeurs  obtenues  indiquent  ainsi  que  la 
superposition  generalement  adoptee  lors  des  relations  structure-activite  de  ces  deux  composes  est 
tres  peu  probable.  A  partir  de  ces  resultats,  les  principes  etablis  du  mode  d’action  des  enkephalines 
au  niveau  du  r&epteur  morphinique  poutraient  se  trouver  considerablement  modifies 
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Conformation  de  molecules  dans  les  membranes: 
RMN  d'enktphalines  et  du  cycloartenol 

A.  Milon3),  T.  Higashijima*3),  T.  Miyazawa^),  J.P.  Kintzinger0),  Y.  Nakatani3),  et 
G.  Ourisson3) 

a)  UA  31 ,  CNRS,  Centre  de  Neurochimie,  67084,  Strasbourg 

b)  Dept.  Biophys.  Biochem.,  Fac.  of  science,  Univ.  of  Tokyo,  Japon 

c)  Laboratoire  de  RMN,  CNRS,  4  rue  Blaise  Pascal,  67008,  Strasbourg 


Nous  presentons  deux  applications  de  la  RMN  a  I'etude  conforma- 
tionelle  de  molecules  interagissant  avec  des  bicouches  lipidiques. 

Le  cycloartenol  est  un  precurseur  biosynthetique  des  phyto¬ 
sterols.  De  nombreux  resultats  in  vitro  et  in  vivo  montrent  que  ce  sterol  a 
des  proprietes  membranaires  indentiques  a  celles  du  cholesterol:  c’est  un 
bon  renforgateur  membranaire,  contrairement  au  lanosterol.  En  vue 
d'expliquer  ces  resultats  il  convient  de  determiner  les  conformations  de  ces 
composes  dans  les  membranes.  Alors  que  le  cholesterol  a  un  squelette 
tetracyclique  rigide,  nous  avons  demontre  par  RMN  et  mecanique  moleculaire 
que  le  cycloartenol  est  une  molecule  flexible  au  niveau  du  cycle  C.  Apres 
attribution  complete  des  spectres  et  13C,  des  experiences  de  COSY  phasee 
a  deux  temperatures  ont  montre  I'existence  d'un  tquilibre  conformationel  en 
solution  0).  Les  resultats  de  RMN,  en  accord  avec  la  modelisation,  montrent 
que  la  conformation  determinee  k  I'etat  solide  (cristallographie)  n'est  pas  la 
seule  existante. 

Les  enkephalines  sent  de  petits  neuropeptides  hydrosolubles  en 
echange  rapide  avec  les  membranes,  pour  lesquels  la  technique  du  NOE 
transfere  en  presence  de  lipides  perdeuteries  a  ete  appliquee  avec  succes. 
Nous  avons  determine  les  conformations  membranaires  de  trois  analogues 
actifs  etcelle  de  un  analogue  inactif  d'enkephalines  (Tyr-LAIa-Gly-Phe-Leu). 

Les  trois  analogues  actifs  ont  la  meme  conformation  membranaire 
constitute  d'un  p-turn  de  type  II'  autour  des  quatre  derniers  residus. 

L'analogue  inactif  a  une  conformation  differente  ©lie  aussi  entierement 
caracterisee  et  qui  permet  de  comprendre  son  inactivite.  Ceci  met  en 
evidence  une  correlation  entre  activite  et  conformation  membranaire.  Ces 
resultats  sont  confrontes  aux  resultats  obtenus  par  dynamique  moleculaire  et 
par  mecanique  moleculaire  (logiciel  SYBYL,  TRIPOS  ass.). 

(1)  A.  Milon,  Y.  Nakatani,  J.P.  Kintzinger,  et  G.  Ourisson,  Helv.  Chim.  Acta,  72, 
1-13  (1989). 
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Etude  th6orique  de  la  structure  du  complexe 
Glutathion  -  Eau  oxyg6n4e 


J.Berg&s.J.Caillet 

Dynamiqiie  des  Interactions  Moleculaires 
Universi^e  Pierre  et  Marie  Curie 

Z.Abedinzadeh.M. Gardes 

Laboratoire  de  Chimie  Physique 
Universite  Rene  Descartes 


Le  glutathion,  en  tant  que  reducteur  intraceilulaire, 
joue  un  role  important  de  protection  de  la  cellule  contre 
differents  agents  toxiques  tels  que  certains  radicaux  libres, 
des  oxygenes  reactifs  ( en  particulier  ceux  des  peroxydes  et 
plus  precisement  du  peroxyde  d'hydrogene)  et  des  composes 
d'origine  endogene  ou  exogene.  II  est  connu  que  la  reaction 
d'oxydation  du  glutathion  par  l'eau  oxygenee  est  catalysee, 
in  vivo,  par  une  enzyme,  la  glutathion  peroxydase.  II  a  ete 
montre  recemment  ( A . Abedinzadeh  et  al.  Can.  J.  Chem .  (1989) 
sous  presse)  que  meme  en  1 ' absence  d' enzyme  le  glutathion 
reagit  rapidement  avec  l'eau  oxygenee  pour  donner  un  complexe 
peroxyde.  La  stabilite  d'un  tel  complexe  est  liee  a  la  for¬ 
mation  de  liaisons  hydrogene  entre  le  glutathion  et  l'eau 
oxygenee  (hypotheses  emises  par  les  auteurs  precedents). 

Nous  avons  done  entrepris  une  etude  theorique  sur  la 
structure  de  ce  complexe  peroxyde,  a  1'etat  isole  ou  en 
presence  de  molecules  d'eau,  apres  avoir  effectue  une  analyse 
conformationnel le  du  glutathion  dans  ses  differentes  formes, 
neutre  et  ioniques.  Les  calculs  ont  ete  faits  par  la  methode 
SIBFA  qui  permet  d'obtenir  simultanement  les  variations 
d' energies  intra  et  intermoleculaires  a  partir  de  formules 
semi -empi riques  ajustees  sur  des  calculs  ab  initio. 
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Etude  Theorique  du  Transfer!  de  Proton  Intramol6culaire 

dans  la  Glycine 

J .  Langlet,  J.  Caillet,  E.  Evleth,  et  E . Kassab 

lynamique  des  Interactions  Ms  ecu  la',  res 
'Jniversite  Pierre  et  Marie  Curie 

Nous  presen tons  cans  ce  travail  les  premiers  resultats  de  calc u Is 
ab  initio  du  chemin  de  transfers  de  proton  m  t  ramo  1  ecu  1  a  i  r  e  dans  la 
q 1  vc  i  p.5  en  t  r  e  ses  deux  formes  n.eut  r  e  ,  N'H,  CK,  COOK  ,  e  t  zwi  t  te  r  i  on  i  cue  , 
NHyCHjCOC  ,  aussi  bien  pour  la  molecule  isolee  que  pour  la  molecule 


solvates 

par  plus: 

sur  s 

me  1  ecu',  es 

d ' eau  . 

Cans  le  cas  d 

e  la 

mo 

1  ecu  1  e 

i solSe ,  1 

e  ca 1 cu 1 

ab  i 

n i t i o ,  au 

niveau 

6-3 1G*  .  mentre 

que 

la 

forme 

es 

t  plus  3t 

able 

que  la  forme  zwi: 

ter i on i que  (de 

1  '  o  r  d 

re 

de  2  5 

Kcal/mol),  et  que  les  minima  qui  leur  correspondent  sur  la  surface  de 
transfer:  sont  separes  par  une  tres  f a ibl e  barriere  ; point 
d ' inflexion) . 

Ncs  calculs  dans  la  meme  base  montrent  que  .ors  de  1  '  hydra  t  a  t  i  on 
pa  r  ceux  itic  .ecu.es  c'eau,  la  forme  neutre  est  tou  jours  la  cl  us  stable 
mais  que  1'ecart  er.erqetique  entre  les  deux  formes  dimmue,  et  passe 
de  25  a  15  Kcal/mol 

1 '  influence  de  la  base  a  ete  egaiement  examinee,  et  notamment 
1'effet  de  fonctions  diffuses  et  de  polarisation  sur  la  stability  des 
deux  formes  neutre  et  zwi  t  ter  iomque  de  la  molecule,  cette  derniere 
s'etant  reveiee  tres  sensible  au  type  de  la  base  utilises. 
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ENVIRONMENT  EFFECT  ON  THE  DOUBLE  PROTON  TRANSFER  IN  CARBOXYLIC  ACID 

DIMERS 

H.  Chojnacki,  M.J.  pyka,  W.A.  Sokaiski 
Inst  i  tute  of  Organic  and  Physical  CLea.  istry  ,  I-- 

Wyb.  Wyspiansk lego  27,  50-070  Wroclaw 
•  POLAND 

s:or.lesf  r*v'ict  ior.s.  ” r. a  double  oro’.on  transfer'  (DPT)  within  the 
s  ubmo  i  ecu  1  tr  :n’ ion  is  assumed  to  h~  .--sponsible  for  some 
biological  processes ,  however.  the  most  conclusive  results  are  still 
iv.i  liable  tor  s  ;  m p  1  e  » o '"lei  s y stems  . 

Quantum  chemical  studies  ot  tautomeric  equilibria  in  isolated 
:  irboxvl ic  acid  dimers  point  out  on  apparent  discrepancy  between  theory 
and  experiment.  In  the  case  of  benzoic  acid  dimer  the  calculated 
-..a  i  vat  ion  energy  tor  DPT  was  estimated  tc  be  5  <y.  al/ac  ie  (or  more)  [1] 
wnereas  its  experimental  value  was  found  to  be  of  the  order  of  1  —  2 
kra  i/niole  [  2— 3  ]  . 

In  this  situation  some  reasons  of  the  above-mentioned  divergence 
must  be  taken  into  account.  One  of  them  is  the  crystal  field  in  the 
solid  state  which  in  the  ordered  lattice  reduces  the  barrier  height  by 
"a.  15  %  [4],  According  to  our  earlier  quantum  chemical  calculations 

[5],  it  is  concluded  that  the  optimal  field  for  the  proton 
rearrangement  may  be  realized  in  unordered  crystal  phases.  The  optimal 
tautomeric  forms  for  the  molecules  around  the  central  dimer  have  been 
predicted  by  expressing  the  environmental  effects  in  terms  of 
i n t ermo 1  ecu  1 ar  interactions  within  the  electrostatic  approximation  [6]. 

•Some  other  sources  of  the  observed  discrepancy  between  experiment 
md  relevant  quantum  ch-micul  calculations  far  fne  proton  transfer  are 
discussed  herein  15,7-3',. 
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EFFETS  STRUCTURELS  ST  DE  CHARGE  SUR  LE  PHTHALOCYANINE  DIMERE 

Francisco  Torrens,  Enrique  Sanchez-Perez  et  Jose  Sanchez-Marin 

Departaraent  de  Quimica-Fisica,  Facultac  de  Quimica,  Universicat 
de  Valencia.  Dr.  Moliner-50,  46100-Bur jassot ,  Valencia,  Spain 


L’agregation  esc  an  phenomena  bien  connu  dans  la  chimie  de 
la  phchalocyanine .  Des  interactions  peuvent  avoir  lieu  entre 
anneaux  de  phthalocyanine  adjacents  aucant  en  phase  organique 
comme  aqueuse.  Ces  associations  resultent  en  1 ' accouplement  de 
deux,  ou  plus,  unices  de  phchalocyanine. 

Un  etude  theorique  de  1 ’ interact  ion  intermoleculaire  entre 
deux  molecules  de  phthalocyanine  a  ete  realise  en  emploiyant  un 
potentiel  de  paires  atome-acome  propose  par  S.  Fraga  [S.  Fraga, 
J.  Comput.  Chem. ,  3  (1982)  329;  S.  Fraga,  Comput.  Phys. 
Commun.  29  (1983)  3511.  Plusieurs  geometries  de  preuve  sont 
optimisees  avec  un  methode  de  metrique  variable. 

Quelques  calculs  preliminai res  montrent  une  structure 
face-a-face  comme  la  conformation  la  plus  stable  pour  le  dimere. 
Dans  cet  arrangement,  une  molecule  peut  etre  tournee  dans  le  plan 
par  rapport  a  1 'autre. 

C'est  bien  connu  que,  en  phase  solide,  une  petite  deviation 
des  positions  moleculaires  relatives  dans  1 ' arrangement  du 
cristal  cause  les  plus  grandes  consequences  sur  les  propcietes 
physiques  du  cristal.  Par  consequent,  nous  avons  etudie 
1* importance  de  quelques  parametres  geometriques  (distance 
interplanaire,  translation  horizontale  des  anneaux,  ...)  de  meme 
que  l'etat  d'oxydation  de  toutes  les  deux  molecules  sur 
1 ‘orientation  relative  des  molecules  individuelles  dans  le 
dimere . 

Aussi,  on  etudie  d'aucres  aspects  du  dimere.  L '  interact  ion 
entre  deux  molecules  apolaires  vient  determinee,  f ondamentalment , 
par  les  moments  quadrupolaires.  On  compare  plusieurs  methodes 
numeriques  pour  le  calcui  des  moments  quadrupolaires.  L' effect 
de  ces  moments  sur  la  geometrie  du  dimere  est  discute. 

Quelques  conclusions  preliminaires  sont:  1)  On  recommande 
1' usage  de  charges  nettes  "ab  initio"  pour  eviter  problemes  de 
t ransf er ibi 1 i tee  des  parametres  d '  interact  ion .  2)  On  recommande 
la  renormalization  des  charges  nettes  pour  reproduire  bons 
moments  quadrupolaires.  3)  L 1 interpolat ion  de  polarizabilitees 
atomiques  resulte  un  procede  approprie  doncs  les  charges  nettes 
sont  modifiees. 
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Th*  inn  solvation  in  water-hexamethy lphospho- 
r  ictr  1  amide  solution:  MD  simulation. 


I .  I .  Vo  l s ma rt ,  M .  6 .  K l s 0  1  o v ,  V.P.  Puh o  vs fc  i ,  V .  M .  Ke s s  i e r 
Institute  of  Non-Pqueous  Solutions  Chemistry 

of  the  USSR  Hcadewy  uf  Sr.  is  no  os ,  Ivanovo,  USSP 


'he  MD  si  my  lotion  technique  has  been  applied 


to  studo  the  solvation  of  the  Na  ion  in  binary 


so  1  ve n  t  water  -he  xame  thy  1  pho  sphor  ic tr  ia.ni  l  de  (HMPT .) . 
C>'»  leu  lot  ions  were  carried  out  on  a  periodic  cub  ic 
box  containing;  one  ion,  one  HMPT  arid  205  rigid  SPC 
water  ruo lecules .  Program  pack;. age  MOOVS  designed  by 
a uth o r s  wa  s  u  s  e d  f o  r  s l mu  1 a 1 1  o  n . 

Five  5-ps  MD  runs  in  NVT-ensemb le  w i th  different 
initial  positions  of  the  ion  were  performed.  Both 
time -aver aged  and  time -dependent  properties  of  wa¬ 
ter  rno  1  e c u  1  e s  we r  e  c a  1  c  u  1  a  ted  i ft  th e  r  e g l o n s  o f 

*  £.  1  -j  +m  4  +U.  r.  r. 1  ►-  4  ».  ►,  _r.  .*»  1  v 

b  i  •  •».  >.b  >J  ‘.‘J  ‘Jib  j.‘'J  i.Jl  ■  1 II  U  M  I  I  j-’ I  •3  I 


site s  of  HMPT,  and  in  the  first,  second  and  third 
s o lvati o n  shell s  o f  th e  l o r i  s e p a r a  tel y . 

Various  radial  and  orientational  distributions, 
v e l o c i t y  a u t o -  a n d  c r oss-co r r e 1  a tion  f u ncti o ns, 
s  tr-  uc  tu  r  e  o  f  h  y  d  r  o  ge  n  -b  o  n  d  e  d  n  e  two  r  k  ,  an  d  th  e  r  mo  - 
dynamic  properties  0*  the  simulated  system  wee* 
studied.  Specific  results  of  mutual  influence  of 
sodium  and  HMPT  on  water  structure  and  tf»e  way  of 
ionic  diffusion  in  solution  were  observed. 
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Monte  Carlo  simulation  of  the 
Acetate-Guanidinium  ion  pair  in  water 


Stephane  30UD0N',  Georges  WIPFF 

Instifut  de  Chimie 

1,  Rue  Blaise  Pascal 

67000  STRASBOURG 

France 


Monte  Carlo  simulations  have  been  carried  out  on  the 
guanidinium/  acetate  ion  pair  in  water  at  25°C  and  1  atm.  The  free 
energy  profile  of  the  separation  of  the  two  ions  has  been  calculated 
using  the  statistical  perturbation  theory.  The  determined  otential  of 
mean  force  between  the  two  ions  set  in  the  configuration  giving  the 
best  interactions  (double-bound  C2V  configuration)  shows  three 
minima.  There  is  no  clear  preference  for  the  intimate  pair  (C  ...  C  = 
3.3A)  The  intermediate  state  (C  ...  C  =  6.3A)  is  a  result  of  cooperative 
binding  between  the  ions  and  two  water  molecules  forming  a  double 
hydrogen-bond  bridge.  The  two  extreme  states  (C  ...  C  =3  .3A  and  9.6A 
respectively)  are  separated  by  a  large  energetical  barrier.  The  computed 
results  provide  insights  in  the  effects  of  the  solvent  on  an  important 
biochemical  system  and  anchoring  site  of  proteins. 
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HIGH  TEMPERATURE  ANNEALED  MOLECULAR  DYNAMICS 
SIMULATIONS  AS  A  TOOL  FOR  CONFORMATIONAL 

SAMPLING: 

A  TEST  CASE  ON  THE  222  BICYCLIC  CRYPT  AND. 


P.  AUFFINGER,  G.WIPFF' 

Laboratoire  de  RMN  et  de  Modelisation  Moleculaire 
L'A  422  du  CNRS 
Institut  de  Chimie 
4  rue  Blaise-Pascal 
67000  STRASBOURG 


ABSTRACT: 


We  have  performed  methodological  studies  using  "High 
Temperature  Annealed  Molecular  Dynamics  Simulations"  (HTAMDS)  on  a 
bicyclic  cryptand,  with  the  following  protocol  :  100  ps  of  molecular 

dynamics  at  1000°K,  followed  by  optimisation  of  the  structures  saved 
every  0.2  ps,  relaxation  of  these  structures  during  20  ps  of  molecular 
dynamics  at  300°K  and,  finally,  optimisation.  Four  sets  of  500  low 
energy  conformers  of  the  free  222  cryptand  have  been  produced 

starting  either  from  the  free  cryptand  or  from  the  M+/222  cryptate  with 

different  representation  of  M+. 

The  analysis  of  these  four  sets  allows  assessment  of  the  ability  of 
this  HTAMDS  technique  :  (i)  to  interconvert  experimentally  known 
conformers  starting  from  one  of  them,  (ii)  to  locate  the  energy  minima, 
fiii)  to  generate  new  conformers  of  low  energy,  and  (iv)  to  account  for 
the  average  structure  observed  on  the  NMR  time  scale.  In  view  of  the 
ionophoric  behavior  of  the  cryptand.  structures  are  analyzed  in  terms  of 
the  "in/out"  orientation  of  the  binding  sites. 

It  is  found  that  the  annealed  simulations  on  the  free  molecule, 
although  sampling  largely  the  conformational  space,  do  not  give 

structures  adequate  for  cation  inclusion;  however  they  generate  the 
lowest  energy  structure  known  experimentally  and  other  new  closely 
related  ones.  Inclusion  of  cation  in  the  simulation  (either  as  a  purely 
electrostatic  "driver”,  or  as  a  charged  sphere)  leads  to  conformations 
found  in  several  complexes. 


2  8  C 


KINETIC  MODELLING  OF  HETEROGENEOUS-CATALYZED  REACTIONS  WITH  THE  ANACIN 
SOFTWARE .  APPLICATION  TO  THE  HYDRODENITROGENATION  OF  PHENANTHRIDINE. 

J.  JOFFRE1 ,  P.  GENESTE1,  A.  GUIOA1,  C.  MOREAU1,  and  G.  SZABO2 
1.  ENSCM,  8,  Rue  Ecole  Normale,  34075  Montpellier,  Cedex  France 
2.  TOTAL-B.P. :  27-  76700  Harfleur,  France. 


'he  hydrodeni trcgenat’ on  of  neavy  hydroca'oons  is  reprinted  or  j 
laboratory  scale  by  a  simple  test  consisting  of  the  conversion  of 
a 1 <y lanilines  in  the  presence  of  nitrogenated  polyaromatic  compounds  lire 
quinoline  cr  phenanthridi''e.  Tress  compounds  have  seen  shown  to  inhibit 
tne  conversion  of  alxylam 1 ines  as  due  to  cneir  important  concentration  in 
the  feed  or  to  the  slow  orogress  in  their  own  conversion.  It  is  quite 
obvious  that  the  innibition  results  from  competitive  adsorption  phenomena 
oetween  tne  substituted  anilines  and  tne  neavier  N-compounds  present  in 
the  feed. 

in  order  to  gain  a  better  knowledge  on  the  innibiting  effects  and 
tneir  possiole  quantification,  a  computer  tool  was  required  allowing  to 
interpret  the  experimental  results  and  to  model  these  innibiting  effects. 

A  preprogrammed  xinetic  model  consisting  --  si,<  reactants, 
intermediates  or  reaction  products  Y1.Y2...Y6  was  defined  : 
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in  whicn  Ad)  are  tne  rate  constants  for  the  forward  reactions,  C'i;  are 
those  for  tne  reverse  reactions  and  b(i)  are  tne  constants  of  adsorption 
for  the  different  compounds  (i). 

'he  reactions  are  assumed  to  be  fi-st  order  in  each  compound  {i ) . 
The  system  of  differential  equations  is  integrated  numerically  by  the 
Runge-Kutta  method.  The  constants  A(i),  C(i)  and  o ( i )  are  optimized  with 
tne  simplex  method  to  fit  in  with  cne  experimental  data. 

in  the  case  of  the  hydrodemtrogenation  of  pnenanthridine,  the 
residue  ootained  after  optimization  of  the  forward  rate  constants  only,  is 
relatively  high  ;  this  residue  becomes  smaller  after  optimization  of  both 
forward  and  reverse  rate  constants.  Finally,  the  smallest  residue  is 
obtained  when  the  optimization  also  takes  into  account  the  adsorption 
constants  of  the  different  compounds  (i). 

Thus,  the  appropriate  software  allows  to  model  the  inhibiting 
effects  of  nitrogenated  polyaromatic  compounds  on  the  conversion  of 
substituted  anilines.  This  software  works  at  a  reduced  running  cost  and 
requires  only  widely  used  computer  equipment.  It  allows  to  determine 
automatical ly  twenty  constants  from  one  set  of  experimental  data  and  to 
visualize  the  results  by  comparing  experimental  and  sitmdated  curves. 
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Le  calcul  de  la  correction  Epstein-Nesbet  au  second  ordre  de  la 
densite  electronique  en  tant  ju'app I ication  repartie  dans  le  contexts 
o'un  reseau  in  format!  cue  local  (  "Network  Computing 1 


Jean-Marie  Lecierca 

-aboratoire  de  Dvnamique  des  interactions  Moieculaires 
Unite  oroore  de  Recnerche  A027I  du  CNRS 

Unite  liee  Dar  Convention  a  1'Universite  Pierre  et  ^arie  Curie  (Paris  5) 


Le  Reseau  infcrmatique  Local  de  ta  Physico-Chimie  Moleculaire  de  1'Universite 
Pierre  et  Marie  Curie  (Paris  6)  constitue  un  espace  privilegie  pour  le  developpement 
d'appl ications  reparties  en  Physico-Chimie  Guantique.  II  est  presente  un  tel 
developpement  pour  le  calcul  de  la  correction  Epstein-Nesbet  au  second  ordre  de  la 
oensite  electronique  dans  le  cadre  suivant : 

-  developpement  du  type  Rayleigh-Schrddinger  , 

-  partition  du  type  Eostein-Nesbet ; 

-  correction  limitee  au  2eme  ordre  , 

-  reference  mono-determinantale  a  couches  completes  , 

-  a)  repartition  sur  une  station  de  travail  multi  modules  de  traitement  ,  ^ 

b)  repartition  sur  plusieurs  stations  de  travail. 

Les  avantages  (en  particular  sur  le  plan  du  rapport  performances/prix)  et  les 
limites  d'une  telle  approche,  par  rapport  a  I'uti lisation  plus  classique  de 
"mainframes''  mono  ou  multi-processeurs,  sont  discutes. 
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Molecular  graphics  for  the  Macintosh. 


Jean-Michel  CENSE 

Laboratoire  d'informatique  chimique 
Ecole  Nationale  Superieure  de  Chimie  de  Paris 
1 1.  rue  Pierre  et  Marie  Curie  -  75231  Paris  Cedex  05  FRANCE 


Now  that  Apple  Macintosh  and  high  resolution  laser  printers  are  widely  available  it  appears  interesting  to 
have  an  easy  access  to  specifically  tailored  packages  allowing  creation  and  manipulation  of  chemical  drawings 
for  direct  printing  or  for  insertion  into  word  processors,  desktop  publishers,  presentation  managers  or  graphics 
software. 

MolDraw  and  MolView  programs  allow  the  production  of  a  wide  variety  of  display  modes  for  3-D  molecules 
(Ball  and  stick.  Space  filling,  Dreiding  and  stereo  pairs  of  all  of  these).  These  3-D  molecules  are  obtained  from 
ASCII  tiles  containing  Cartesian  coordinates,  fractional  coordinates  or  internal  coordinates. 

MolDraw  and  MolView  programs  read  texfiles  of  atomic  coordinates  according  to  several  formats,  SYBYL 
PDB  and  a  private  format  for  Cartesian,  fractional  or  internal  coordinates.  These  text  files  are  transfered  from 
host  computer  or  created  on  the  Macintosh. 

Cristallographic  files  include  cell  constants,  fractional  coordinates,  symmetry  operators  and  duplicate 
operators  to  display  several  contiguous  cells. 

Several  parameters,  such  as  covalent  radii,  atomic  radii,  patterns  and  gray  tones  are  adjustable.  Using  a 
Macintosh  ft  and  a  color  monitor  all  models  are  displayed  in  color. 

X,  Y  and  Z  rotations,  translations,  inversions,  scaling,  rotations  around  a  non-ring  bond  are  provided.  An 
axis  defined  by  to  atoms  can  be  moved  to  X,  Y  or  Z  axis.  A  plane  defined  by  three  atoms  can  be  moved  to  XY, 
YZ  or  ZX  plane.  A  cristallographic  molecule  may  be  projected  onto  a  "hid"  plane. 

Information  about  interatomic  distances,  angles  or  torsions  are  available  by  selecting  pertinent  atoms. 

Several  techniques  have  been  tested  for  displaying  depth  relationships  :  depth-cueing,  stereopsis  and 
animation. 

Stereo  views  are  projected  in  either  a  relaxed  mode,  a  mirror  mode  (one  view  is  seen  after  reflexion  on  a 
mirror)  or  a  crossed  mode.  The  stereo  separation  and  the  rotation  between  stereo  views  are  adjustable. 

Animation  is  a  very  effective  technique  for  displaying  depth  relationships.  Successive  views  of  a  3-D  model, 
projected  after  rotation  around  any  axis,  are  stored  in  memory  and  later  transferred  to  the  screen.  While  it  is 
impossible  to  infer  the  sense  of  rotation  from  orthonormal  projections  of  wire-frame  models,  using  depth- 
cueing,  ball  and  stick  or  space  filling  models,  the  sense  of  rotation  can  be  easily  deduced.  The  number  of 
animated  images  and  the  animation  speed  depend  on  the  size  of  the  front  window  and  on  the  size  of  available 
memory. 

Most  models  are  generated  in  less  than  3  seconds  on  screen.  A  black  and  white  shaded  space-filling  image  of 
a  500  atoms  molecule  was  generated  in  30  seconds  on  a  Macintosh  II  at  screen  resolution.  The  same  colored 
image  was  generated  in  15  seconds  only.  Shaded  3-D  models  are  printed  in  less  than  three  minutes  at  maximum 
resolution  (300  ppi). 
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Phase  behaviour  of  cyclohexane  fieri  ved  from  differential  thermal 
analysis  and  computed  hv  molecular  dynamics 

*'v  f! .  Hrheisel,  Theoretisohe  Chenie 

and  a  .  Vilrfl  inner,  Phvsikal ische  Ch.emie  II 

n  i :  hr  -!  n  i  vn  r  s  i  t  p  1  R nc h i :n  ,  ■ir  3  Roc  hi  :r  ,  F  ’’0 

i  fferent  i  al  thermal  analysis  lias  been  nset!  to  study  the  phase  ho- 
tnvioiiPS  of  cy  1  obevene  tn_l'  nnr  deliberated  c\  r  !  o  hp\  an  n  '  r  ^  . 

u  o  t  h  substances  ex  h  i  n  i  t  n  1  as  tic  arr*  non  — rotator  n  t  p  c  n  «  p«  up  1  1  pc 
n  e«  h i rh -nr es si  ire  inducer  obese  transitions, 

A  six-centre  l.ennarr1- Pones  potential  is  nrncnser1  for  eve  1  rhev  ,-re 
anr<  tester!  bv  rolecular  dynamics  ( VP )  calculations.  The  thermrdv- 
nflnic  and  transport  properties  of  experimental  0^1-'  in  liouid  ant 
solid  states  were  reproduced  in  acceptable  approximation.  In  par¬ 
ticular  the  phase  transition  to  the  hi  oh  temperature  rotator  phase 
< althotirh  amorphous)  is  correct.lv  predicted  by  the  calculation. 
Furthermore  the  first  two  rear i en ta t i on  a  1  correlation  functions  of 
three  vectors  coincidinn  with  different  directions  of  the  molecule 
have  been  invest i nated,  from  which  valuable  insinht  into  the 
rotational  behaviour  of  the  occurrinn  linuid  and  solid  phases  is 
derived. 
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ELECTROSTATIC  PROPERTIES  OF  SOLVATED  PROTEINS:  A  MICRO¬ 
SCOPIC  ANALYSIS  BASED  ON  COMPUTER  SIMULATIONS 
Shoshana  J.  Wodak,  Daniel  Van  Belle*  and  Martine  Privost, 

Unite  de  Conformation  de  Macromolecules  Biologiques,  Universite  Libre  de  Bruxelles, 
and  * Plant  Genetic  Systems,  CP160,  P2,  Ave.  P.  Heger,  B-I050  Bruxelles,  BEL¬ 
GIUM 


Molecular  dynamics  simulations  of  a  small  protein  bamase  in  presence  of  ex¬ 
plicit  water  molecules  are  described  and  results  from  a  100  ps  trajectory  of  the 
system  are  analyzed. 

The  deviations  of  the  average  protein  conformation  from  the  starting  crystal 
structure  is  very  satisfactory  (-1  Arms  for  backbone  atoms),  and  the  agreement 
between  computed  and  crystallographic  atomic  fluctuations  is  unusually  good  for 
portions  of  the  protein  that  do  not  participate  in  crystal  contacts.  The  structure 
of  water  around  polar  and  non-poiar  groups  on  the  protein  surface  also  seems 
reasonable  in  that  it  agrees  well  with  previous  observations  made  in  hydrated 
crystals  or  in  simulations  of  small  solvated  molecules. 


The  thrust  of  our  study  concerns  the  use  of  the  last  25  ps  of  the  generated  tra¬ 
jectories  to  obtain  a  detailed  microscopic  description  of  electrostatic  interactions 
in  the  protein.  Contributions  to  these  interactions  from  permanent  protein  di¬ 
poles,  from  orientable  solvent  dipoles,  and  from  electronic  polarizability  are 
evaluated.  Although  the  analysis  must  still  be  qualified  as  preliminary,  a  number 
of  clear  trends  emerge.  The  contribution  of  water  to  local  fields  inside  the  pro¬ 
tein  is  substantial.  It  affects  field  magnitudes  and  field  orientation  respectively 
by  30%  and  31  degrees  on  the  average.  In  comparison,  the  contribution  from 
electronic  polarizability  alone  is  much  lower  with  an  average  of  6%  in  field 
magnitude  and  12  degrees  in  field  orientation.  But  statistical  analysis  of  both 
contributions  shows  that  they  display  an  appreciable  degree  of  inhomogeneity 
throughout  the  protein  matrix  which  suggests  that  their  relative  importance  may 
vary  according  to  the  local  environment. 


Modelling  of  Molecular  Structures  and 
Properties  in  Physical  Chemistry  and 
Biophysics 
Nancy  -  1989 
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NMR  studies  of  proteins  and  protein-DNA  interactions 


J.A.C.Rullmann, 

Department  of  Organic  Chemistry,  University  of  Utrecht, 
Padualaan  8,  3584  CH  Utrecht,  the  Netherlands. 


With  two-dimensional  NMR  spectroscopy  solution  structures  of  small  and  medium 
sized  biomolecular  systems  can  be  determined.  The  nuclear  Overhauser  effects,  or 
NOE's,  in  the  biomolecule  are  translated  into  proton-proton  distance  constraints,  which 
can  be  used  in  Distance  Geometry  algorithms  and  Restrained  Molecular  Dynamics 
simulations  to  determine  the  range  of  possible  structures.  The  accuracy  of  the  distances 
derived  from  NOE's  can  be  improved  by  an  iterative  relaxation  matrix  approach.  The 
procedure  will  be  illustrated  by  the  structure  determination  of  DNA  octamers  and  of  the 
protein  crambin. 

A  major  bottleneck  for  determining  structures  by  NMR  is  the  assignment  of  proton 
resonances.  By  performing  three-dimensional  experiments  the  resolution  can  be 
increased.  Some  results  of  non-selective  three-dimensional  measurements  will  be 
shown. 

Two-dimensional  NOE  studies  of  complexes  of  lac  repressor  headpiece  with  lac 
operator  fragments  are  presented,  including  a  2:1  complex  with  a  22  basepair  operator 
(total  mol.wt.  25.000).  A  large  number  of  protein-DNA  NOE's  has  been  identified.  By  a 
combination  of  molecular  docking  and  Restrained  Molecular  Dynamics  using  these 
NOE's  as  constraints,  models  of  the  lac  headpiece  operator  complex  were  built.  An 
interesting  feature  of  the  structure  of  the  complex  is  that  the  second  ("recognition")  helix 
binds  in  the  major  groove  of  the  DNA  with  an  orientation  that  is  approximately  180  0 
different  from  what  has  been  found  for  othe,  DNA  binding  proteins.  The  model  predicts 
specific  protein-DNA  interactions  which  have  been  found  recently  in  genetic 
experiments  with  the  complete  lac  repressor  as  well. 
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Protein  structure  determination  using  NMR  data  with 
energy  minimisation  and  molecular  dynamics  calculation 

A.  MIKOU,  E.  GUITTET,  V.STOVEN,  M.A.  DELSUC, 
and  J.Y  LALLEMAND. 

Institut  de  Chimie  des  Substances  Naturelles 
CNRS,  91190  Gif  sur  Yvette,  FRANCE. 


Recent  developments  in  two  dimensional  NMR  spectroscopy 
allow  the  structure  determination  of  proteins  in  solution. 

These  structures  are  calculated  by  a  combination  of  NOESY  * 
determined  distances  with  energy  minimisation  and  molecular 
dynamics.  These  distances  are  introduced  as  constraints  during  the 
calculations. 

However,  before  these  methods  can  be  applied,  it  is  necessary  to 
assign  the  various  resonances  from  the  NOE  maps  to  their  parent 
protons. 

The  M.C.D  strategy  (  Main-Chain-Directed  )  (1)  for  assigning 
these  resonances  is  used. 

The  quantitative  evaluation  of  2D  maps  using  the  "  P.A.R.I.S  " 

(  Program  for  Automatic  Recognation  and  Integration  of  2D  NMR 
Signals  )  (2)  provides  information  of  interproton  distances  smaller 
than  5  A. 

This  approach  is  illustrated  with  the  toxin  from  AaH  scorpion 
venom. 


(1)  Englander,  S.W.,  Wand,  A.J.  (1987)  Biochemistry  26,  5953-5958. 

(2)  Stoven,  V.,  Mikou,  A.,  Piveteau,  D.,  Guittet,  E.,  and  Lallemand,  J.Y. 

(1989)  J.Magn.Reson.  in  press. 


*  NOESY:  Nuclear  Overhauser  Effect  Spectroscopy. 
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INFLUENCE  OF  THE  BASE  SEQUENCE  AND  CONFORMATION  ON  THE 
STRUCTURE  OF  THE  POLYNUCLEOTIDE  HYDRATION  SHELL 

Y.I. Poltev,  A.V.Teplukhin,  G.G. Malenkov 

Institute  of  Biological  physics,  USSR  Academy  of  Sciences, 
Pushchino,  Moscow  Region,  USSR 

The  hydration  shells  in  the  minor  and  major  grooves  of 
several  double  helices  with  different  nucleotide  sequences  and 
conformations  have  been  studied.  The  Monte  Carlo  method  (algo¬ 
rithm  of  Metropolis  et  al. )  'was  used  for  simulation  of  systems 
containing  helical  stacks  of  complementary  base  pairs  and  30 
water  molecules  per  each  base  pair  with  unidinensional  periodic 
boundary  conditions  along  the  axis  of  the  helix.  The  repeating 
unit  comprised  6  base  pairs-.  The  mutual  arrangement  of  the  base 
pairs  corresponded  to  the  A  and  3  conformations  of  p-  p(iA). 
poly(dT),  poly(dG) ‘poly(dC) ,  poly(dA-dT) •pcly(dA-dJ)  and 
poly ( dG-dC )• poly (dG-dC) .  The  average  energetic  and  structural 
characteristics  for  the  simulated  systems  at  room  temperature 
and  at  a  temperature  close  to  absolute  zero  were  calculated, 
the  patterns  of  the  hydrogen  bonded  bridges  formed  by  1,2  and 
3  water  molecules  connecting  hydrophilic  centres  of  the  bases, 
are  considered  and  the  probabilities  of  formation  of  3uch 
bridges  have  estimated.  The  hydration  shell  structure  depends 
not  only  on  the  nature  of  the  base  pair,  but  on  the  sequence 
and  conformation  although  the  total  hydration  characteristic- 
(energy  and  its  components,  the  number  of  water-water  and  water- 
-base  H-bcnds)  depend  only  slightly  on  the  latter  factors.  The 
structural  elements  peculiar  to  the  hydration  shells  of  the 
minor  and  major  grooves  of  helical  stacks  of  base  pairs  have  been 
found.  The  probability  of  formation  of  one  water  molecule  bridges 
between  the  hydrophilic  centres  was  found  to  increase  when  the 
stack  transforms  from  the  A  to  3  configuration  (one  water  mole¬ 
cule  bridge  means  that  a  water  molecule  forms  H-bonds  with  two 
or  three  base  atoms).  This  result  is  discussed  in  connection 
with  the  water  economy  concept  proposed  for  the  explanation  of 
the  B  to  A  DNA  transition.  The  difference  in  the  structure  of  the 
hydration  shells  near  NHg-groups  of  A:T  and  G:C  pairs  was  found 
and  discussed. 
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SIMULATION  OF  PROTEINS: 
STRUCTURES,  DYNAMICS,  &  THERMODYNAMICS 

•Martin  Karplus 
Department  of  Chemistry 
Harvard  University 


Proteins  are  one  of  the  essentlaL  constituents  of  living 
systems  so  that  an  understanding  of  their  properties  Is 
necessary  for  a  fundamental  approach  to  bloLoglcal  problems. 

It  will  be  demonstrated  that  molecular  dynamics  simulations 
provide  a  method  for  obtaining  Information  concerning  protein 
structure,  dynamics,  and  thermodynaml cs •  Of  particular 
Interest  are  the  Internal  motions  which  are  sizeable  In  spite  of 
the  close  packing  of  the  atoms  In  the  native  structure;  such 
fluctuations  pLay  an  Important  role  In  ligand  binding.  The 
utlLlty  of  simulations  for  structure  refinement  and  for 
determining  the  properties  of  modified  proteins  will  be  described. 
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Base  Sequence  Effects  and  Transitions  in  DNA 


Richard  Lavery 

Xnstitut  de  Biologie  Physico-Chimique 
13,  rue  Pierre  et  Marie  Curie 
Paris  75005,  France 


Experimental  evidence  on  the  biological  role  of  fine 
structure  within  DNA  is  rapidly  accumulating.  In  principle, 
theoretical  modeling  should  be  able  to  help  in  deciphering  this 
new  "code",  however  the  size  and  complexity  of  the  molecular 
systems  involved  has  hindered  progress  in  this  area. 

Recent  developments  in  our  laboratory  have  led  to  two  new 
methodologies  which  should  improve  the  situation.  Firstly,  we 
have  developed  an  energy  minimisation  procedure,  specifically 
oriented  to  the  treatment  of  nucleic  acids,  which  directly  uses 
helicoidal  parameters  as  variables.  We  are  thus  able  to  describe 
DNA  oligomers  with  ten  times  fewer  variables  than  are  necessary 
in  classical  molecular  mechanics.  At  the  same  time  we  are  able  to 
study  much  more  easily  the  energy  dependence  of  structural 
deformations.  Secondly,  we  have  formulated  an  algorithm  for 
rigorously  describing  the  conformation  of  irregular  nucleic  acid 
oligomers  and,  in  particular,  their  curvature.  Application  of 
these  techniques  to  studying  the  influence  of  base  sequence  on 
the  fine  structure,  the  flexibility  and  the  conformational 
transitions  of  DNA  will  be  presented. 
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ETUDE  CONFORMATIONNELLE  PAR  DYNAMIOUE  MQLECULAERE 
DE  QUELQUES  CnHIBITEURS  DES  COLLAGENASES  BACTERIENNES. 

V.  DIVE  1.  B.  GELQUIN  !.  D.  PERAHIA  2 


1  Service  de  Biochimie,  Departement  de  Biologie.  CEN'-Sac!av,  Gif-sur-Yvette 

2  Laboratoire  d'Enzymoiogie  Physico-Chimique  et  Moleculaire,  Bat  430. 

University  Pans-Sud,  9I4Q5  Orsay  Ceaex. 


Les  collagenases  bacteriennes  sont  des  metailoproteines  a  zinc  capa'oles 
de  diver  avec  une  grande  specificite  le  collagene  natif  et  des  substrats 
de  synthese.  Une  famille  d’inhibiteurs  de  I’enzvme.  de  formule  generals 
HS-CH2-CH2-CO-PRO-X  (X=  L-aminoacide),  a  ete  obtenue.  Les  constables 
d'inhibition  mesurees  varient  fortement  avec  la  nature  du  residu  X. 

Les  derives  N-methyles  de  X  presentent  une  activite  plus  grande  par 

rapport  au  compose  parent.  Ils  sont  caracterises  par  une  reduction  de 
l'espace  conformationnel  par  rapport  au  compose  non  methyle.  Deux  de 
ces  composes  (X  =  Ala  et  X  =  N-methly-Ala)  ont  ete  choisis  en  vue  de 
1'etude  par  dynamique  moleculaire  de  la  restriction  de  l'espace 

conformationnel . 

Les  calculs  de  dynamique  moleculaire  realises  a  600K  a  l'aide  du  logiciel 
CHARMM,  ont  permis  de  mettre  en  evidence  1'existence  de  3  families  de 
conformations  pour  le  compose  normal  et  de  2  families  pour  ie  compose 
methyle.  Les  valeurs  des  barrieres  de  potentiei  et  des  minima  d'energie 
pour  ces  different  conformeres  ont  ete  calcules  en  geometrie  relaxee. 
La  comparaison  .ec  les  resultats  experimentaux  de  la  RMN  montre  la 
necessite  de  la  prise  en  compte  dans  ces  calculs  du  facteur  entropique  et 

des  effets  de  solvam. 

Dans  ce  but.  plusieurs  sous-programmes  ont  ete  developpes  et  integres 
au  logiciel  CHARMM.  11s  cermettent  de  caiculer  la  difference  d'energie 
libre  entre  differents  conformeres  a  partir  de  dvnamiques  moleculaires 
a  300  K. 

La  prise  en  compte  du  facteur  entropique  donne  un  bon  accord  avec  les 

resultats  experimentaux  de  la  RMN. 

Ce  travail  est  poursuivi  actuellement  pour  tenir  compte  des  effets  de 
solvant. 


EVALUATION  OF  ELECTROSTATIC  PROPERTIES  AT  ENZYME  ACTIVE  SITES. 


G.Dive,  J . Lamotte-Brasseur  and  D.Dehareng 

Service  de  Microbiologie ,  University  de  LiSge, 
Sart  Tilman  , Belgium. 


Understanding  at  the  atonic  level  how  proteins  and  their 
ligands  interact  is  one  of  the  most  crucial  challenges  in  the 
molecular  interaction  studies.  Such  an  interaction  is  charac¬ 
terized  by  an  energy  E(int)  that  may  be  evaluated  at  different 
levels  of  approximation  depending  on  the  size  of  the  system  and 
the  distance  between  the  interacting  partners. 

At  large  distances,  the  most  important  component  of  E(int) 
is  the  electrostatic  interaction  energy  E(es).  This  work:  attempt 
to  calculate  E(es)  from  valuable  expressions  of  the  electrostati 
potential  (E.P.)  at  the  quantum  chemistry  level. 

In  the  first  part,  modelized  active  sites  of  alpha  chymo- 
trypsin  and  of  native  and  mutated  subtilisins  have  been  inves¬ 
tigated  in  view  to  point  out  the  electrostatic  features  involved 
in  the  beginning  of  the  interaction.  The  level  of  approximation 
in  the  derivation  of  the  E.P.  factors  (density  matrix  elements, 
monoelectronic  integrals),  the  influence  of  the  number  of 
residues  in  the  model  (from  91  to  277  atoms),  as  well  as  the 
influence  of  the  solvent  molecules  (from  0  to  25  water  mole¬ 
cules)  are  discussed.  The  most  striking  result  is  the  relation 
between  the  direction  of  the  high  dipole  moment  value  and  an 
impressive  E.P.  negative  well  located  on  one  side  of  the  cleft 
in  and  around  the  nucleophilic  serine  CProteins,  J.L.B.  et  al, 
submitted] . 

The  second  part  of  this  work  proposes  a  numerical  procedure 
to  derive  E(es)  with  a  high  performance  ratio  (accuracy  level)/ 
(CPU  time).  This  procedure  consists  in  a  3D  numerical  integra¬ 
tion  of  the  product  of  the  E.P.  of  one  partner  and  the  charge 
density  of  the  other.  The  approximations  involved  in  this 
calculation  refer  to  the  way  the  electrostatic  potential  and 
the  charge  density  are  determined.  Numerical  ab  initio  SCF 
ST0-3G(Wls)  as  well  as  CNDO  type  calculations  are  oresented 
versus  analytical  SCF  ST0-3G(Wls)  ones.  These  clearly  show  that 
the  numerical  procedure  at  the  semiempirical  level  is  quite 
satisfactory  as  far  as  the  accuracy  is  concerned,  and  has  the 
great  advantage  to  be  far  less  time  consuming  than  the  ab  initio 
analytical  way  CTh.Chim. Acta,  D.D.  et  al,  accepted]. 


Theoretical  Prediction  of  Base  Sequence  Effects  in  DNA  and  RNA 

Brigitte  HARTMANN 

Institut  de  Biologie  Physico-Chimique 
13,  rue  Pierre  et  Marie  Curie 
Paris  75005,  France 


The  role  of  the  base  sequence  in  determining  the  structure 
of  3-DNA  is  already  -.veil  established.  But  little  is  known  about 
base  sequence  effects  for  other  conformation  such  as  Z  or  A 
allomcrphs . 

The  use  of  a  new  molecular  modeling  methodology,  JUMNA, 
based  on  internal  and  helicoidal  variables,  has  allowed  us  to 
detect  clear  effects  of  base  sequence  on  the  conformation  of 
DNA  and  RNA  double  helix. 

The  results  we  shall  present  for  oligomers  with 
dinucleotide  symetrv  imposed  show  that,  especially  within  the 
Z-DNA  or  RNA  family,  considerable  polymorphism  can  exist. 
Certain  sequences  can  adopt  more  than  one  conformation.  In 
these  cases,  the  helicoidal  characteristics  are  very  different 
although  the  energies  are  very  similar.  We  note  that,  overall, 
large  changes  in  base  pair  positioning  require  only  small 
coupled  changes  in  backbone  dihedral  angles. 

Comparison  between  B-  and  Z-DNA  or  A-  and  Z-RNA 
conformations  show  that  the  order  of  B-Z  transition  enthalpies 
generally  correlates  with  the  number  of  syn-purine/anti- 
pyrimidine  errors  (i.e.  syn-pyrimidine  and/or  anti-purine). 

Finally,  we  use  our  molecular  modelling  for  the 
interpretation  of  the  highly  variable  experimental  chemical 
reactivities  of  a  natural  Z-DNA,  in  terms  of  conformational 
discontinuities . 
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Modeling  DNA  Curvature  as  a  Function  of  3ase  Sequence 


Marc  Poncin 

Insfitut  de  3iolcgia  Physico-Chirique 
*3  2TM3  ?i.22rr,3  ■2‘t.  Ci lr’is 

Paris  ”5005,  France 


It  is  now  well  known  that  certain  sequences  of  IN' A  are 
naturally  curved  and  that  this  curvature  can  play  an  important 
role  in  the  biological  functioning  of  DMA.  However,  today  there 
is  still  no  clear  structural  explanation  for  these  effects. 

We  have  consequently  undertaken  a  theoretical  modeling 
study  in  an  attempt  to  better  understand  hew  the  base  sequence 
can  influence  the  fine  structure  of  DNA.  The  results  presented 
will  show  some  applications  of  the  programs  Jumna  and  Curves 
which  have  been  specifically  designed  for  the  construction, 
energy  minimisation,  manipulation  and  analysis  of  nucleic  acid 
oligomers.  The  use  of  helicoidal  parameters  and  internal 
coordinates  during  energy  minimisation  is  particularly  helpful 
for  investigating  base  sequence  effects  and  DMA  flexibility  in  a 
controlled  way.  The  rigorous  analysis  provided  by  the  Curves 
algorithm  enables  us  to  study  both  natural  curvature  and  to 
investigate  the  anisotropy  of  imposed  DMA  bending. 


38  A 


MINIMISATION  ET  DYNAMIQUE  MOLECULAIRE  DE  L'ADN-Z 
MODI  FIE  PAR  L'ACETYLAMINOFLUORENE 


V.  FP.ITSCH  et  E .  WESTKOF 

Laboratoire  de  Cristallographie  Biologiaue 
Ir.stitut  de  Biologie  Moleculaire  et  Cellulaire, 
Centre  National  de  la  Recherche  Scientif ique, 
15  rue  R.  Descartes,  F-67034  Strasbourg. 


Le  programme  de  mecanique  moleculaire  AMBER,  version 
3.0,  a  ete  utilise  pour  etudier  les  conformations  adoptees 
par  la  guanosine  modifiee  sur  la  position  CS  par  le 
carcincgene  acetylaminof luorene  (AAF)  .  A  cet  effet,  nous 
avons  introduit  dans  le  programme  AMBER  une  routine 
permettant  d'utiliser  une  "constante  dielect r ique 11  prenant 
des  valeurs  entre  1  et  80  ;  celles-ci  varient  de  fagon 
sigmoidale  en  fonction  de  la  distance,  suivant  la  forme 
suggeree  par  Lavery,  et  al.  (J.  3iom.  Struct.  and 
Dyn . , 3 , 989-1  Oil  (1986  )). 

Nous  avons  effectue  une  minimisation  systematique  de 
la  guanosine  modifiee  par  l'AAF  en  fonction  de  l'angle  de 
torsion  autour  de  la  liaison  glycosidique  et  de  celui  autour 
de  la  junction  entre  la  guanosine  et  l'AAF.  Les 
conformations  optimales  ainsi  obtenues,  compatibles  avec  la 
forme  Z  de  1'ADN,  ont  ete  inserees  dans  un  hexamere-Z 
d'CGCGCG>2  et  a  nouveau  minimisees.  Les  structures 
minimisees  ont  ensuite  ete  soumises  a  un  calcul  de  dynamique 
moleculaire  partielle  sur  la  sequence  5 ' -C-GAAF-C-3 '  .  Tous 
ces  calculs  ont  ete  effectues  avec  une  "constante 
dielectrique"  egale  a  r,  a  4r,  ou  "sigmoidale". 

Une  des  conformations  preferees  aussi  bien  en 
minimisation  qu'en  dynamique  est  representee  dans  la  figure 
s te reoscopique  ci-dessous. 
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RECONNAISSANCE  DES  BASES  DE  L'ADN  PAR  LIAISON  HYDROGENE 
Effet  de  la  longueur  d'une  chaine  hydrocarbonee 

par  C.  COULOMBEAU(1>  et  N.  GRESH<2) 

Universite  J.  Fourier,  LEDSS  VI,  B.P.  53F,  38041  Grenoble 

Cedex,  France. 

Institut  de  Biologie  Physico-Chimique ,  13,  rue  Pierre  et 

Marie  Curie,  75005  Paris,  France. 


OBJET.  L’ne  serie  de  derives  du  methoxy-2  Chloro-6  amino  9 
acridine  pcrtant  une  chaine  aliphatique  de  longueur  variable 
terminee  par  un  groupe  carboxamide  libre,  interagit  de  maniere 
preferentielle  avec  un  ADN  riche  en  paires  (G-C).  Cette 
interaction  varie  avec  la  longueur  de  la  chaine.  On  estime 
qu'elle  est  optimum  lorsque  le  groupe  terminal  peut  former  dans 
le  petit  sillon  deux  liaisons  hydrogene  avec  le  groupe  NH2  et 
le  N3  d'une  guanine  adjacente  au  site  d' intercalation. 

RESULTAT  DES  CALCULS .  Le  calcul  par  SIBFA  de  ces  derives  de 
l' acridine  intercales  dans  un  hexamere  rigide  forme  de  paires 
de  bases  (G-C)  montre  qua  la  formation  d'un  complexe  bidente 
n'est  pas  indifferent  a  la  sequence  :  ils  se  forment  dans  une 
sequence  (GCCGGC)2  mais  pas  dans  (CGCGCG)2. 

Pour  ces  complexes  bidentes  l'energie  d' interaction 
stabilisante  est  la  plus  importante  pour  le  compose  comportant 
5  methylenes,  ce  qui  est  conforme  a  1 ' experience . 

De  meme  la  difference  entre  l'energie  totale  du  complexe 
et  l'energie  du  ligand  seul  etendu  indique  que  la  stabilisation 
est  maximum  pour  l'isomere  a  5  methylenes.  La  comparaison  des 
energies  des  complexes  pour  des  sequences  d'hexameres 
( GCCGGC) 2 ,  (ATTAAT)2  et  (TATATA) 2  preserve  le  caractere 
favorable  pour  le  meme  compose. 

On  peut  noter  cependant  que  ces  sequences  particul ieres 
sont  insuf f isantes  pour  rendre  totalement  compte  des  faits 
experimentaux . 
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MOLECULAR  MODELLING  OF  TWO  DOUBLE  STRANDED  DNA  FRAGMENTS 
TGACGTCA and  ACTGCAGT 


J.M.  PIRIOU1,  O.  MAUFFRET1,  J.  ARMIER2.  S.  FERMANDJIAN1.  M.  LE  BRET1 

1  Laboratoire  de  Biochimie-Enzymologie  CN'RS  URA  158,  INSERM  U140 
Institut  Gustave- Roussy-  39,  rue  Camille  Desmoulins 
94805  Villejuif  Cedex  FRANCE 

"  Institut  de  Recherches  Scientifiques  sur  le  Cancer 
7,  rue  G.  Moquet  -  94800  Villejuif  FRANCE 


Footprinting  experiments  showed  that  the  antitumour  drug  N-methyl  9-hydrcxy 
ellipticine  (NMHE)  preferrentially  binds  to  CG  sites  in  double  stranded  DNA.  and  led  to  the 
determination  of  a  consensus  sequence  ACGT  (B.  RENE  and  co-workers). 

To  elucidate  this  selectivity,  we  have  carried  out  the  *H  NMR  study  and  the 
molecular  modelling  of  the  binding  of  NMHE  to  a  DNA  octamer  containing  the  consensus 
sequence  -  TGaCGTCA  -,  and  its  inverse  sequence  ACTGCAGT.  The  two  fragments  have 
first  been  studied  without  the  drug. 

On  the  basis  of  these  NMR  data,  we  undertook  the  molecular  modelling  of  the 
oc  tamers. 

From  the  coupling  constants  Jl’2\  J1’2”,  J2'3’,  J2"3’,  measured  from  COSY 
experiments,  we  have  considered  each  sugar  ring  as  the  result  of  an  equilibrium  between  the 
C2’-endo  and  the  C3’-endo  conformations,  which  can  be  represented  by  the  percentage  of 
C2’-endo  form  (Rinkel  and  Altona  1987). 

Thus  for  both  molecules  a  set  of  64  =  26  starting  structures  were  generated, 
considering  for  each  sugar  ring  (except  those  at  the  extremities)  the  C2'-endo  and  C3'-endo 
conformations.  Those  models  were  refined,  without  constraining  the  geometry  of  the 
deoxyribose.  We  selected  the  best  ten  ones  in  terms  of  energy  (with  less  than  5  kcal 
difference  between  one  another).  We  calculated  the  coupling  constants  of  our  models  and 
confronted  them  to  the  experimental  constants.  We  considered  each  nucleotide  independently. 
For  some  of  them,  the  calculated  couplings  were  in  good  agreement  with  the  experimental 
ones  (pure  C2’-endo).  For  the  other  ones  we  had  to  consider  the  equilibrium  C2’-endo  -  C3’- 
endo,  i.e.  taking  the  C2’-endo  form  of  one  sugar,  we  tried  to  find  another  structure  where 
the  same  sugar  exists  in  the  C3’-endo  conformation,  and  a  fraction  p  such  as,  for  the  4 
couplings  : 


P"1C2'endcj  +  1 1  ~PVc3'«ndp 

the  structures  tor  whi 


=  J 


exp 

it  was 


not  possible  to  find  such  a  fit  for  all 


We  ruled  out  the  structures  for  which 
nucleotides.  Thus  we  kept  3  models. 

Then  interproton  distances  were  used  for  the  final  determination  of  the  structures. 

The  work  concerning  the  binding  of  the  drug  to  the  fragments  is  currently  being 


done. 
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NUCLEOTIDE  MISPAIRS  STABILIZED  3Y  NATSR  BRIDGES. 

MODELING  OF  STRUCTURE  AND  ROLE  IN  TEMPLATE  BIOSYNTHESIS 

7. I.poltev,  3. V. Steinberg 

Institute  of  Biological  Physios  USSR  Academy  of  Sciences, 
pushchino,  Moscow  Region,  142292,  USSR 


To  elucidate  the  mechanisms  of  wobble  in  codon-anticodon 
interaction  and  of  errors  in  template  biosynthesis  of  nucleic 
acids,  the  possibility  of  insertion  of  water  molecules  into 
wobble  nucleotide  pairs  was  considered.  At  the  first  step  the 


calculations  of  intermclecular  inte¬ 

raction  energy  for  systems  containing  two  bases  and  one  or  two 
water  molecules  were  carried  out  by  the  method  of  atom-atom 
potential  functions.  There  exist  energy  minima  for  each  base 
pair,  corresponding  to  a  single  N-H...C  or  N-R...N  H-bond  bet¬ 
ween  the  bases  and  H-bonding  of  the  water  molecules  with  both 
bases.  In  the  other  minima,  where  two  bases  are  connected  via 
two  H-bonds,  water  bridges  additionally  stabilise  such  pairs. 

The  position  of  the  bases  and  water  molecules  of  seme  calculated 
complexes  correspond  to  those  of  nucleotide  pairs  in  some  double 
-helical  oligonucleotides,  as  known  from  X-ray  and  NMR  data. 

At  the  second  step  we  have  attempted  to  insert  these  pairs 
into  the  third  position  of  the  codon-anticodon  complex  -under 
the  assumption  that  the  rearrangement  of  the  sugar-phesphate 


backbone,  that  i3  necessary  for  such  a  pair  formation,  occurs 
only  in  the  anticodon,  but  three  nucleotides  of  codon  during 
codon-ant  icodor.  recognition  have  the  fixed  conformation  of  A-RNA 
Consequently,  along  with  the  standard  pairs  the  following  pairs 
become  probable:  UG,  'JI,  UU,  CU,  GU  and  AI,  GU  and  AI  being 
weaker  than  the  other  ones.  Ail  other  pairs  appear  to  be  steri- 
cally  impossible.  Atom-atom  calculations  of  the  system  (codon- 
-anticodon  loop  -  -water  molecule)  established  the  energy  minima 
for  every  permitted  wobble  pair  and  revealed  that  the  conforma¬ 
tional  rearrangement  of  the  anticodon  loop  which  accomodates 
the  latter  for  each  of  these  pairs  is  fairly  small  and  so  quite 
possible.  All  the  conclusions  are  in  good  agreement  with  the 
experimental  data. 
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STUDY  ON  THE  SEQUENCE  SELECTIVITY  OF  THE  BI-FUNCTIONAL  INTER CALATOR 
- KTTRCALISTUM - 


R.C.MAROUN.  M.DELEPIERRE  AND  B.P.  ROQUES 

Laboratoire  ass  c onformaticns  et  des  interactions  de  modeles  bio'.cgiques  et 

pharmacologiques 

et 

Unite  de  Pharmacochimie  Moleculaire 
(UA49S  CNRS  et  U266  INSERMi 

4.  av  de  i'Observatoire.  75006  Pans 
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ABSTRACT 


The  structures  and  binding  energetics  of  the  complexes  formed  between  the  rigid 
linking  chain  DNA  bi-functional  intercalator  ditercalinium  (1)  and  a  series  of  model 
tetradecxynucleotibe  duplexes  that  follow  always  an  alternating  punne-pyrimidine 
sequence  have  been  investigated  by  means  of  JUMNA  (2).  a  molecular  mechanics 
approach.  The  ditercalinium  complexes  all  contain  one  internal  excluded  site  such 
that  two  base  pairs  are  sandwiched  between  the  chromophores  intercalated  at  the 
terminal  sites. 

A  comparative  energy  analysis  of  all  the  complexes  has  permitted  us  to  order  them 
according  to  their  base  pair  preference.  The  energy  minimisation  calculations  using 
the  JUMNA  methodology  point  to  an  increased  stabilisation  of  the  minor  groove 
complex  (i.e..  the  complex  in  which  the  piperidine  rings  of  the  linking  chains  of  the 
drug  occupy  the  minor  groove  side)  over  the  corresponding  major  groove  complex  in 
the  case  of  d(ApTpApT)2.  All  other  complexes  are  major  groove  complexes.  These 
results  put  forward  the  many-sidedness  of  the  modes  of  intercalation  of 
ditercalinium.  We  have,  in  addition,  characterised  quantitatively  the  local 
deformations  induced  by  the  drug  on  the  DNA  fragment. 

The  prediction  in  the  case  of  the  ATAT  sequence  is  in  agreement  with  the  theory  of  the 
negative  molecular  electrostatic  potential  of  the  DNA  minor  groove  proposed  by- 
Pullman.  Lavery  and  coworkers. 

(The  mode  of  intercalation  of  ditercalinium  with  the  ApTpApT  sequence  is  presently- 
being  investigated  in  this  same  laboratory  by  1 2H  and  3 Ip  NMR  techniques.) 


(1)  Pelaprat,  D.,  Delbarre,  A,  Le  Guen.  I..  Roques,  B.P.  and  Le  Pecq.  J.B.  (1980)  J.  Med. 
Chem  22-  1336-1343. 

(2)  Lavery.  R  (1988)  in  'DNA  Bending  and  Curvature,"  W.K.01son.  M.  Sundaralingam. 
M.H.Sarma  and  RH.Sarma,  Eds.  Adenine  Press.  New  York.  NY.  191-201. 
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DNA  Bis- Intercalation  in  the  7H  pyrido  [4,3-cl  carb azole 
family  :  A  comparative  1H  and  31p  NMR  study. 


Muriel  DELEPIERRB\  Rachid  MAROUN *,  Catherine  AGLHET*.  Christiane  GARBAY 
Jean  IGOLEX+,  Bernard  ROgUES *. 


*  Laboratoire  de  chimie  organique  ,  CNRS  UA  498, 1XSERM  L’266.  PARIS  75006. 
+  Unite  de  chimie  organique,  CNRS  UA  4S9.  Institut  PASTEUR- 


Ditercallnlum  (2,2’-u4.4'-biplperidineI- 1.  l'-diyldi  2.1  ethane  dlyl)bis(10- 
methoxv  -  7H  pyrido  [4.3C]  carbazoliumltetramethane  sulfonate)  (NSC  366241)  a  DNA 
bis-intercalating  compound,  is  a  potent  antitumoral  rigid  dimer.  Previous  studies  (1) 
have  shown  that  a  reduced  flexibility  of  the  linking  chain  of  such  dimer  is  essential 
for  its  biological  activity.  In  order  to  understand  at  the  molecular  level  the 
mechanism  of  action  and  the  structure-activity  relationship  of  these  series  of  DNA 
lntercalators.  new  dimers  with  additional  methylene  groups  between  the  two 
piperidine  rings  have  been  synthesized  (2).  Addition  of  one  methylene  group  in  the 
chain  preserved  the  activity,  whereas  addition  of  two  methylene  groups  reduced  the 
cytoxlcity  which  finally  disappears  when  three  methylene  groups  are  inserted. 
Therefore,  the  study  of  the  interaction  of  dimers  bearing  no  (202)  two  (222)  and  three 
(232)  methylene  groups  with  the  self  complementary  hexanucleotide  d(CGATCC)2  have 
been  Investigated  by  1H  and  31p  NMR 

The  results  indicate  that  all  dimers  bislntercalate  into  the  DNA  and  exhibit  a 
neat  preference  for  pyrimidine-purine  sequences.  The  intermolecular  NOE  effects 
between  the  dimers  and  the  nucleotide  lead  to  the  conclusion  that  the  three  dimers 
intercalate  with  their  rigid  bls-ethyl  blspiperidine  chain  fitting  the  major  groove  of 
the  helix.  Inter-residue  NOE  effects  at  the  DNA  level,  as  well  as  induced  shifts  are 
discussed  in  relation  to  both  the  conformational  changes  induced  on  DNA  upon 
intercalation  and  to  the  different  activity  of  the  dimers. 
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BIS -INTERCALATION  OF  DITERCALINIUM  AND  "FLEXIBLE" 
ANALOGUE  IN  THE  OCTANUCLEOTIDE  d(TpTpCpGpCpGpApA)?.: 

A  COMPARATIVE  STUDY  BY  NMR  AND  MOLECULAR  MODELLING 


Joel  POTHIER  'A  Muriel  DELEPIERRE®  ,  Marie-Christine 
BARSI®.  Christiane  GARBAY®  .  Jean  IGOLEN*  .  Marc  LE  BRET1-  & 
Bernard  ROQUES®  . 


-  Laboratoire  de  Physico-Chimie  Macromoleculaire  ,  URA  158  CN'RS. 
Instltut  Gustave  Roussy.  rue  Camille  Desmoulins.  94S00  VILLEJUIF. 

8  Departement  de  chimie  Organique  U266  1NSERM.  UA  498  CN’RS,  UER  des 
Sciences  Pharmaceutiques  et  Biologiques.  4  Avenue  de  l'Observatoire.  75006  PARIS 
France. 

;  Unite  de  chimie  Organique.  Institut  Pasteur.28  rue  du  Docteur  ROUX 
75015  PARIS  France 

•  To  whom  correspondance  should  be  addressed 


7H-Pyrido  [4,3  ]c  carbazole  dimers  of  the  ditercalinium 
family  are  DNA  bis-intercalators  that  display  high  DNA  affinity  and 
antitumor  properties.  The  latter  depends  on  substitution  on  the 
pyridocarbazole  rings  and/or  "flexibility''  of  the  linking  chain.  When 
the  two  piperidine  rings  of  ditercalinium  (  figure  I  a  )  are  replaced 
by  six  methylene  groups  (  as  shown  in  figure  I  b),  the  drug 
cytotoxicity  desappears. 

The  interaction  of  Ditercalinium  and  its  "flexible"  analogue 
with  both  the  tetranucleotide  d(CpGpCpG)2  and  the  octanucleotide 
d(TpTpCpGpCpGpApA)2  has  been  therefore  investigated  by  'H  and 
31P  NMR.  The  octanucleotide  was  chosen  as  it  has  only  one  site 
available  for  bis-intercalation.  and  is  long  enough  to  probe  long 
distance  perturbations. 

The  results  indicate  that  all  dimers  bis-intercalate  into  the 
DNA  via  the  major  groove  with  similar  geometry  but  with 
substantial  differences  at  the  DNA  level,  particularly  in  the  sugar- 
phosphate  backbone. 

In  addition,  molecular  modelling  calculations  are  being 
carried  out  on  these  complexes  and  compared  to  the  models 
obtained  from  NMR  data. 
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DYNAMIQUE  MOLECULAIRE  POUR  LA  SIMULATION 
DE  REACTIONS  PONT  LA  CINETIQUE  EST 
LIMITEE  PAR  LA  DIFFUSION 


W.  Dong,  F.  Boros  et  J.C.  Andre 


GRAPP-QCPR-UA  323  du  CNRS,  ENSIC-iNPL,  1,  me  Grandville  BP  451.  5-1001  NANCY  CEDEX 


Les  constantes  de  vitesse  apparente  d'un  grand  nombre  de  reactions  en  phase  liquids  sor.: 
affectees  par  les  phenomenes  de  transport  moleculaire  De  nombreux  modeles.  touiours  Sondes 
sur  l'hypothese  d'un  solvent  considere  comme  un  continuum,  ont  ete  proposes  pour  expliquer  de 
tels  phenomenes.  Ils  presentent  tous  des  imperfections,  plus  ou  moins  flagrantes  suivant  le  type 
de  svsteme  moleculaire  considere. 

^experimentation  directe  est.  bien  sur,  la  seule  possibilite  d'etablir  'ear  vaiidite.  Mais  eile 
n  est  pas  toujours  realisable,  a  cause,  soit  d’une  precision  insuffisante,  soil  du  choix  d  un 
svsteme  moleculaire  adapte.  De  ce  point  de  vue.  la  dynamique  moleculaire  fournit  des 
"experiences  theoriques"  et  de  preeieux  renseigr.ements  sur  la  vaiidite  des  concepts,  mcme  s'il 
s'agit  d'une  methode  un  peu  lourde  et  couteuse.  ne  pouvant  suppleer  aux  modules  "continus" 
lorsqu'on  veut  interpreter  des  resultats  experimentaux. 

Un  calcul  de  dynamique  mole'culaire  a  ete  encrepris  pour  etudier  ces  reactions  dont  la 
cinetique  est  limitee  par  la  diffusion  moleculaire.  Toutes  les  molecules  du  milieu  considere 
(solvant  et  molecules  reactives)  sont  supposees  spheriques,  de  taille  identique,  la  reaction  se 
produisant  des  le  premier  contact  (hypotheses  de  SMOLUCKOWSKI).  Nos  premiers  resultats 
montrent  l’existence  d'effets  transitoires  non  prevus  par  les  modeles  continus.  permettant 
d'expliquer  en  partie,  par  exemple,  les  ecans  observes,  dans  les  graphes  de  STERN  VOLNLER. 
lorsqu'on  excite  de  faqon  continue  (c'est  a  dire  en  integrant  tout  le  declin)  ou  pulsee  (c'est  a  dire 
en  essayant  de  resoudre  temporellement  le  declin). 

De  telles  simulations  peuvent  maintenant  eae  entreprises  dans  le  but  de  decrire  des  svstemes 
plus  complexes  (tailles  differentes  pour  le  solvant  et  les  molecules  reactives,  reactivate  non 
uniforme,  potentiels  iocaux,  etc...). 
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Structure  and  Dynami ca  of  Water  and  Ionic  Solution  near 
Blcmembrane  Surfaces  from  Molecular  Dynamics  Simulations 

yi.  Schler.krich ,  K.  Nlcklas,  Ph.  3opp  and  J.  3rickmanr. 
Instltut  fur  Physikallsche  Chamia 
Technische  Hochschula  Darmstadt ,  Patarsenstr ■  20 
D-6100  Darmstadt,  F.  R.  G. 

The  interface  between  a  biomembrane  surface  and  liquid  water  {or 
ionic  solution)  is  studied  with  molecular  dynamics  simulation 
technique.  Therein  the  membrane  surface  is  modeled  by  CCC“ 
groups  and  a  TIP4P  model  is  used  for  water  molecules .  Strong 
ordering  of  the  water  molecules  close  to  the  membrane  is  found. 
Ions  influence  this  structure  as  well  as  that  of  the  surface. 
The  ion  transport  from  the  bulk  phase  towards  the  membrane 
surface  is  studied. 
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TRANSPORT  DANS  LES  MEMBRANES  BIOLOGIQUES 


MODELISATION  ET  EXPERIENCES 


E.  Garcin1,  F.  Baros1,  J.C.  Andre1,  D.  Daveloose-,  J.  Viret-, 
M.L.  Viriot1  et  M.  Donner3 


1)  GRAPP-DCPR-L'A  323  5u  CNRS,  ENSIC-INPL.  i  rue  Grandvilie,  3?  451. 54C01  Nancy  Cedex  France 

2)  CRSSA,  Laboracoire  de  3iophysique,  BP  37. 38702  La Tronche  Cedex  France 

3)  U  234  INSERM.  Plateau  de  Srabols,  CO  10.  5451 !  Vandceuvre  les  Nancy  Cedex  France 


L  crude  des  membranes  biologiques  er  plus  parricuiieremenr  de  la  dvnamique  membranaire 
fair  appel  a  differentes  techniques  physiques  d'investigarion,  qui  permettent  de  caracteriser  soit  un 
parametre  d'ordre,  soit  les  proprietes  de  transport  en  terms  de  coefficients  de  diffusion  laterals. 

Pour  des  systemes  modeles  simples  (soit  des  phospholipides  purs,  soit  des  melanges 
consumes  de  phospholipides  additionnes  de  proteines  membranaires  purifiees)  nous  avons,  a  la 
t'ois,  realise  des  experiences  et  une  modelisation  qui  mettent  en  evidence  un  role  d’obstacle  de  la 
part  des  proteines  sur  les  proprietes  de  transport  dans  les  membranes. 

L'expenmentadon  repose  sur  1'emplci  de  trots  techniques  complementaires.  la  fluorescence 
resolue  dans  le  temps  (inhibition  de  la  fluorescence  du  pyrene),  la  RPE  (nitroxydes  :  5  et  16  NS) 
et  la  calorimetrie  differenrielle. 

La  modelisation,  qui  ttent  compte  de  concepts  de  physique  des  liquides 
pseudo-bidimensionnels.  permet  de  simuler  les  proprietes  de  transport  au  sein  des  membranes  et 
de  prevoir  l'effet  des  proteines. 


A  SIMPLE  MODEL  OF  NUCLEIC  BASES  ASSOCIATION  :  THE  CYCLIC  DIMER  OF  THl 

'  T^PYETUORE - 

VIBRATIONAL  STUDY 


J.  FAVROT,  Th.  BOUISSOU,  J.F.  BRAZIER,  U.A.  au  CNRS  454,  Laboratoire  Synthese,  Structure  et 
Reactivate  de  Molecules  PhosDhorees 


S.  CASTILLO ,  Laboratoire  de  Recherche  sur  VEnergie 

A.  ZWICK,  U.A.  au  CNRS  74,  Laboratoire  de  Physique  des  Sondes 

Umversne  Paul  Sabatier,  1 1 8  route  de  Narbcnne.  3 !  062  Toulouse  cede* ,  France 


Vibrational  spectroscooy,  specially  ultraviolet  resonance  Raman 
spectroscopy,  is  often  used  to  study  nucleic  bases  association  or  tauto- 
merlsm,  nucleic  scids  structural  chsn^os  or  nuclsic  scids  intersection  with 
proteins,  metal  ions  or  drugs  [Tj.  For  such  studies,  it  woula  be  of  interest 
to  correlate  the  frequency  shifts  or  the  intensity  variations  of  the  lines  to 
the  association  state  of  the  nucleic  bases. 

As  simple  model  system  we  have  chosen  the  equilibrium  between  toe 
2-pyridone  monomer  (M)  and  the  dimer  (D)  of  this  molecule  : 
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By  infrared  spectroscopy,  we  have  found  the  concentration  ranqe 
giving  nearly  the  same  intensity  of  monomer  and  dimer  mq.q  bands 

directly  involved  in  hydrogen  bond  association  [2]  namely  ~  2. 1C'1 2  M  in 
CH3CN  or  CD3CN  and  ~  1 0"3 4  M  in  CHC13  or  CDC13  solutions. 

In; the  present  work  we  try  to  assign  the  secondary  effects  of  hydro- 
gen  ucnding  on  heterocyclic  ring  leading  to  distinct  frequencies  for  tne 
monomer  (free  2-pyridone)  and  for  the  dimer  (associated  2-pyridone). 

previous  work  on  pyridinium,  pyrazinium  or  pyrimidmium  salts  (31  or 
on  2-pynmiaone  and  2-pyrimidone  salts  (4]  have  shown  that  some  rinq 

freauencies  Sr  P  td  t-ha  hnnWind  nf  - X 


or 
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icies  are  sensitive  to  the  bondina  of  an  hydroqen  atom  to  nitroaen 
Nn  ).  With  respect  to  the  soectra'o;'  tne  solid, 'for  2-cvr;dcn-  sole- 

bnmfl  h  o  rvM  c  -a  r\  r\  _ _ i  _  .  _  •  j. _ ■  •  - 


t'ons,  some  pands  appear  as  doublets.  The  relative  intensities  of  each  cor., 
ponent  of  these  doublets  vary  with  concentration.  Then,  the  spectral 
changes  are  related  to  the  evolution  of  the  2  M^  D  eGuilicnum.  Infrared 
spectra  of  2-pyridone  solutions  (2. 1 0~ 1  M  in  CD3CN)  show  such  doublets  at 

L47ki’  '33R~'367,  731-724  and  518-510  cm-1.  The  components  of 
these  doublets  are  assigned  to  monomer  and  out  of  phase  motion  of  the 
dimer,  we  use  also  Raman  spectroscopy  in  order  to  identify  the  in  phase 
motion  of  the  centrosymmetric  dimer. 


1 .  Proceedings  of  the  eleventh  International  Conference  on  Raman  Spectroscopy,  5-9  September 
1 988,  London,  England.  Ed.  by  R.J.H.  CLARK  and  D.A.  LONG,  John  Wiley  and  Sons,  1983 

2.  S.  CASTILLO,  J.  FAVROT ,  Th.  BOUISSOU  et  A.  ZWICK,  Communication  aux  Journees  de  Chimie 
Physique,  1989,  Paris 

3.  R.  FOGLIZZO,  These  de  Doctorat  es  Sciences  Physiques,  1970,  Paris 

4.  E.  PICQUENARD,  These  de  Doctorat  es  Sciences  Physiques,  1 982,  Peris 


i  thj  i  i  •  on  of  Pol /ionic  Interactions  by  Prooressive  Aooroxim&t ions 


Heiko  H.  Ohlenbusch 

Til '.'ERG  I  "PE  LOUIS  PASTEUR  DE  STRASBOURG,  "ACULTE  DE  MEDECINE,  INST  I  TUT  DE 
PKYsTCUE  SIOLOGIQUE .  4,  rue  K i rsch 1 eoer ,  r-S70Sv  STRASBOURG  Cedex 

.emeu  ter  generated  ionic  site  cornel  at  i  on  functions  applied  to  polyionic 


:u  »;  were  anal /zed  in  order  to  obtain  into 
c r oba t- ■  1  i  t y  of  polv ionic  interactions.  It  i 
pe‘", i  ts  predictions  about  preferred  mutual 
sno  could  thus  be  used  as  ou i de  for  as 


mat  ■  on  about  the 
s  shown  that  this 
orientations  of 
time. tine  the  most 


peace  ;>;■  of  approach  of  interacting  charged  oartic'es  as  well  as 
e  : :  -  c  r  -  a  t  -  r.  for  ionic  comp  1  exes  ,  which  corresponds  to  that  of  the 
west  e'e:*.-p  static  free  energy.  ’he  method  o  *  progressive 

or  ox : ma  t  one  can  be  applied  step-wise  to  e 1 ec trostat i c  interactions 
-  then  for  greater  detail  to  the  analysis  of  the  effects  of  dielectric 
encmena.  electric  pc  lari  cat  ion  or  other  short-range  me!  ecu's.'  forces, 
ain  .  r.  3  step -wise  manner  .  Ado  I  i  cat  ion  of  the  analysis  to  highly 
arced  mcdei  systems  has  shown  that  it  is  possible  to  select  geometric 
rre’j*  on  s  between  alternative  poly  ionic  conformations  bv  neolectino 
=t  5*  the  sner t-nange  molecular  forces.  The  results  obtained  sugoest, 
at  the  sf-ects  due  to  solvent  interactions,  counterions  and  electric 
pc  :  ar  :  cat ;  on  are  dominated  by  long-range  electrostatic  interactions  as 
T tr  as  over  all  conformational  adaptation  is  concerned,  when  the 

me! ecu !*», taking  oar t  in  the  interactions,  are  highly  charged  species, 
a  ‘-re  5 f  site  correlation  function  to  be  used  to  estimate  the 

rength  of  the  molecular  interactions  depends  quite  naturally  on  the 


strength, 
t.-pe  of 


to  estimate  the 


of  the  molecular  interactions  depends  quite  naturally  on  the 
interacting  sites  investigated:  van  der  Uaals  interactions  are 


coun ter  1  on; 


a  =  •.  i  iT.a  tec  rom  functions  ot  the  type  f  ( l -'R^o)  ,  interactions  in  the 
presence  cf  counterions  from  f ( l/(Ri j*e*p<k*Ri j > > } ,  the  electrostatic 
■  nferact  ions  frem  the  strength  of  the  electrostatic  field,  i.e. 
e  '•  !  /  •  R  i  .1  *P :  ;  n  etc.,  where  R  and  Ri.j  are  the  interatomic  distances  and 
k>  the.  Oeb  ■  e-Huecke  1  parameter. 

•  h  Fig.  1  p  w;ts  of  l/( R*R)  are  shown  for  a  moving  pel /ionic  probe  and  a 
eo  do  I » . on  in  a  coordinate  system.  The  fly-by  correlation  curves  show 
alternating  maxima  of  repulsions  and  attractions  and  permit  the  pairwise 
alignment  of  the  molecules  corresponding  to  their  interaction  maximum. 
Imoe  the  geometric  information  content  derivable  from  such  fly-by 
correlation  curves  does  not  depend  very  much  on  the  value  of  the 
sxpor.ent  o*  P ,  ‘he  structural  deductions  that  may  be  drawn  from  such 
analyses  b:  not  depend  very  much  on  the  exact  t-.-pe  of  correlation 
•unction  0‘csen.  Th  i  s  alignment  determines,  in  this  part  cular  example, 
tre  Position  of  the  diad  axis  of  symmetry  of  the  polvionic  dimer, 

T1"?  me“-.cb  analysis  has  been  applied  to  complex,  h • ghl  /  charged 

r  acr  orr.pl  ao.  a=.  '.ba  u,;,r  instance  histones,  the  basic  nuclear  proteins. 
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1  s  s  e  n  s  c  h  .a  f  t  e  n 
F inure  1 


Fly-by  correlation  curve  of  an 
artificial  polyion,  polyaianine 
<AI a) 50  with  blocked  ends  and 
G1  u 1 9 , 26 , 33 . 40 , Ar  g22 , 2? , 3 6 , 43 . 
designed  to  mimic  the  central 
segment  of  histone  J. 
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INHIBITOR  AND  SUBSTRATE  BINDING  AND  CATALYSIS  IN  D1HYDROFO- 
LATE  REDUCTASE  -  A  THEORETICAL  APPROACH 


T.  Solmajer 

Boris  Kidric  Institute  of  Chemistry 
Hajdrihova  19 
61000  Ljubljana, 

Slovenia 

Yugoslavia 


A  fundamental  problem  in  chemistry  and  biochemistry  of  proteins  is  understanding  the 
role  of  a  single  amino  acid  residue  in  determining  the  molecular  properties  of  the  protein. 
Site-specific  directed  mutagenesis  is  a  very  powerful  approach  which  yields  experimental 
data  such  as  changes  in  catalytic  properties  for  mutant  species  of  enzymes  compared  with 
their  native  analogues.  Still  more  understanding  can  be  obtained  of  an  enzyme  system  if 
the  X-rav  structures  of  the  parent  and  mutant  enzymes  are  solved.  E.Coli  dihydrofolate 
reductase  (wild  type,  Ser  27  and  Asn  27  mutants)  is  such  a  system  and  experimental  data 
from  ( I)  -  the  X-ray  structures,  binding  constants  of  methotrexate  inhibitor  (MTX)  and 
catalytic  activity  for  substrate  dihydrofolate  -  were  used  as  a  basis  for  our  theoretical  study. 


A  novel  methodology  consisting  of  a  combination  of  ab  initio  split  valence  basis  set  cal¬ 
culations  and  the  thermodynamic  integration  method  implemented  with  the  molecular 
dynamics  was  used  to  determine  the  binding  constants  and  activation  energetics  of  the  in¬ 
hibitor  and  substrate,  respectively.  The  results  obtained  support  the  view  that  the  electros¬ 
tatic  forces  play  a  major  role  in  enzyme  catalyses. 

References: 
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Relationship  between  molecular  properties  of  some 
Photosystem  XI_  herbicide  and  their  activity  on  tnazme- 
resistant  and  susceptible  chloroplasts . 

CREYZZT  3.  v'LQUXM  1  Z  M  D:'CR’-T3XT.  Service 
3iophysique.  Dept  de  3iologie.  OZA/IIIRA  Saclay.  31191  Sif-su 
Yvette . 
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methyl  chiopyri.midine  herbicides  inhibit  the  phucrsy 
electron  flow  by  competing  with  a  pi  astoquir.cr.e  cole:', 
binding  at  a  c or, rr.cn  site  or.  the  secondary  electron  acce 
?  bio  t  o  s  v  s  t  sni  "ZZ. ,  ^ orrr.ed  bv  Z *  -'v- ^ ~ 

The  conformations  and  electronic  properties  of  seve 
these  compounds,  representative  cf  different  str 

families  active  on  the  same  site,  were  determined 
molecular  mechanics  i  SHARMM ,  3Y3Y1,;  or  quantum  semi-em 

:  MNDO ,  AMD  methods.  Dihedral  angles,  rotational 
barriers,  dipole  moment  and  ionisation  potential 
calculated.  Spectroscopic  measurements  of  the  in  vitro 
bonding  and  charge  transfer  capacities  of  these  compound 
also  tarried  tut  ar.d  compared  with,  the  cal  related  cha: 
densities.  Results  obtained  by  AMI  were  found  to  be  the  m 
consistent  with  experimental  data  and  this  method  was  used 
study  the  hydrogen  bonding  properties  of  the  carbonyl  cf  ur; 
or  the  ring  nitrogens  of  trianmes  with 
donors . 

The  conformational 
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ETUDE  STRUCTURALE  COMPARATIVE  D' HYDROXY  ET  DE  METHOXY  FLAVQNES. 
CCMPARAISCN  ENTRE  LE5  RE5ULTAT5  RX  ET  DIVER5E5  METHGDES  DE 
MODEL I SAT  I CN  MOLECULA I  RE . 


j.C.  wALuE-,  P.  BCNZrAoSI*,  -  •  M  .  CAXDCU 

^ 3  r  3 1 ^ 0  ^r'v^ocr^iiTii.s  1  is  i  i  Un ivsrsits 

3  Ai x-Marse 1 1 1  a  3, 

*  iu"  at  Cantos  as  Caicui  de  1  jnivsrsits  au  Mans) 


_0S  *  1 avcncicss  sent  uns  c 1 a555  is  Cw^pcsss  f rsquemmen ~ 

'anccntrSs  dans  las  olantas.  Css  aamoosas,  avac  Zaurs  analogues 
syntnstiquss  sent  a io 1 og iquemen t  actifs.  Mous  nous  semmes 
intaresses  au*  a  ~  *  a  1 3  canrcruationnels  das  interactions  antra 
Za  cycle  pnenyie  at  la  systame  fietsrocyci  isus  cnromcne. 

i„ss  atuoes  RX  de  t-ois  ■‘iavones  net.ncxvlees  ;  la  3-nethoxy,  la 
-metnexy  at  la  2  ,  o  -d  me  tnoxy  )  mettent  an  evidence  une 

-rsicn  cu  cycle  phenyie  variant  de  queiques  eagres  pour  la  2  - 
ne  tnoxy  a  70,7°  pour  la  2  ,  i'-di.T  .-thoxy  flavors. 

Civerses  netnodes  da  mode  1  isation  no  1  ecu  1  ai re  ,  de  type 
Cuantique  i MNDO ,  AMI)  at  de  necanique  no  1  ecu  1  airs  (MMl! 
nontrent  qua  la  metnode  AM  1  ce  M.J.5.  Dewar  donna  las  resultats 
las  cl  us  prcc.nes  de  cau  .  da  1  atat  c-istallm.  _s 
genera  1 isation  aux  f  Uvones  pol ymStha*y less  centenues  dans  1 es 
Citrus  pernettra  d'envisager  des  correlations  structurs- 
ac  1 1  v  i  te  . 
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Structure  determination  from  NMR  using  a  relaxation  matrix  approach  and 
stereospecific  assignment  procedures;  application  to  the  solution  structure  of  crambin. 

J.A.C.Rullmann,  M.M.G.Koning,  R.M.J.N.Lamericns,  R.Boelens  and  R.Kaptein, 
Deoartment  of  Organic  Chemistry,  University  of  Utrecht, 

Paduaiaan  3,  358-1  CH  Utrecnt,  the  Netherlands. 

An  important  prcoiem  in  structure  cetermination  using  2D  NMR  data  is  the  accurate 
evaluation  of  distance  constraints  from  observed  proton-proton  NOE  cross  peaks. 
Ideally,  the  distance  can  be  derived  from  the  initial  buildup  rate  of  NOE  intensity  as  a 
function  of  the  mixing  time  used  in  the  NMR  experiment.  Since  at  small  mixing  times  the 
intensity  is  also  small,  the  initial  rate  cannot  be  measured  directly.  At  longer  mixing 
times,  however,  magnetisation  transfer  from  or  to  ether  neighbouring  protons  (the  so 
called  spin  diffusion)  has  to  be  taken  into  account.  We  describe  the  Iterative  Relaxation 
Matrix  Approach  (or  IRMA  procedure),  wnich  handles  spin  diffusion  quantitatively  by 
solving  the  Bloch  equation  for  the  NOE  magnetisation  transfer  simultaneously  for  all 
spins  in  the  system.  After  combining  experimental  and  theoretical  NOE  intensities  (the 
latter  being  derived  from  a  model  structure),  the  relaxation  matnx  desenbing  the  kinetics 
of  the  proolem  is  constructed.  From  this  the  spin-diffusion  corrected  constraint  aistances 
can  be  derived  directly.  Using  restrained  Molecular  Dynamics  an  improved  structure  is 
calculated,  which  in  turn  is  used  to  derive  an  improved  relaxation  matrix.  The  whole 
procedure  can  be  repeated  until  all  experimental  NOE’s  are  explained  satisfactorily. 
Theoretical  NOE’s,  i.e.  cross-peak  intensities  predicted  on  the  basis  of  a  structural 
model,  can  also  be  used  for  stereospecific  assignment  of  proton  pairs  which  are 
otherwise  difficult  to  discriminate,  e.g.  31/62  protons  and  terminal  methyl  groups  in  Val 
and  Leu.  It  will  be  shown  how  these  assignment  procedures  can  be  implemented,  and 
how  they  can  be  combined  with  structure  optimizations  in  which  the  barrier  to 
interchange  of  the  equivalent  protons  has  been  reduced  to  zero.  Strategies  for  the 
ootimizaticn  of  structures  with  Distance  Geometry  and  restrained  Mciecular  Dynamics 
are  discussed,  using  the  solution  structure  of  crambin,  a  protein  with  dg  residues,  as  an 
example. 


COMBINED  USE  OF  NMR  AND  RESTRAINED  COMPUTER  SIMULATIONS  FOR  THE 


MODELING  OF  THE  CYCLIC  ANTIFUNGAL  LIPOPEPTIDE  :  STENDOMYCINE . 


M.Genest,  B.Stawarz,  J.P.Simorre,  D.Genest,  A.Caille  and 


M.PtaJc. 


CBM.CNRS.,  1A  Av.  de  la  Recherche  Scientif ique, 
45071  Orleins  Cede;-:  02 


For  determining  the  three  dimensionnal  structure  of  a 
molecule  in  solution,  the  most  powerful  method  is  based  on 
measurements  of  proton-proton  nOe's.  Because  this  mechanism  is 
strongly  distance-dependent ,  nOe's  can  be  translated  into 
proton  proton  distance  used  as  constraints  in  additional  energy 
term  in  computer  simulations. 

But,  translating  nOe's  in  term  of  distances  is  not  always 
unambigous  due  to  internal  molecular  motions  and  and  spin 
diffusion  phenomena. 

Refinement  structures  based  on  the  superposition  of  the 
experimental  and  calculated  nOe's  intensity  matrices  seems  a 
better  strategy  than  the  fitting  of  the  estimated  and 
calculated  distances  matrices  to  take  account  of  all  the 
protons  proximities. 

This  refinement  method  is  applied  to  the  conformational 
analysis  of  a  cyclic  peptide,  stendomycine ,  an  antifungal 
lipopeptide . 
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MOLECULAR  DYNAMICS  STUDIES  OF  PRISTINAMYCIN  II4  :  FROM  CRYSTAL 
STRUCTURE  TO  CONFORMATIONS  IN  APOLAR  SOLVENT  USING  NMR  DATA. 


E.  SURCOUF,  I.MORIZE,  D.FRECHET,  M.  VUILHORGNE, 

A.MIKOU*,  E.GUITTET*  &  J.Y.LALLEMAND* 

Rhone-Poulenc  Same,  Centre  de  Recherches  de  Vitry,  13  Quai  Jules  Guesde,  BP  14, 
94403  Vitry  sur  Seine,  France. 

*CNRS-ICSN.  91190  Gif  sur  Yvette.  France. 


Pristinamycins  I  and  II  are  natural  antibiotics,  first  isolated  from  Streptomyces 
pristinaespiralis.  These  two  groups  of  molecules  are  chemically  different,  as  PI  is  a  peptidic 
macrolactone  and  PII  an  olefinic  macrolactone,  but  they  biologically  act  in  a  synergic  way.  Each 
group  contains  several  closely  related  molecules;  PIIA  is  the  major  component  of  PII. 


An  hemisynthetic  work  has  been  realized  by  Rhone-Poulenc  Sante  to  find  active 
water-soluble  Pristinamycins.  At  the  same  time,  crystallographic  studies  of  these  two  families 
were  initiated. 


In  the  present  work,  the  conformation  of  PIIA  in  CDCL3  was  determined  using  restrained 
molecular  dynamics  simulations  with  the  help  of  proton/proton  distance  constraints  obtained  by 
ID  and  2D  NOE  spectroscopy  at  400  and  250  MHz. 


The  software  package  INSIGHT/DISCOVER,  the  programs  MANOSK  &  Chem-X  were 
used.  The  structure  thus  obtained  was  compared  to  the  X-Ray  crystal  structure  previously 
published  (F.Durant  &  ah.Cryst.Struct.Comm.,  1974,3,503-510).’ 
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Etude  eonf ormationnelle  de  la  tete  polaire  de 
1 ' Amphot her i ci ne  B  &  I'6tat 
isol£  ou  en  presence  de  molecules  d'eau. 


J.BergSs,J.Caillet>J.Langlet 

Dynamique  des  Interactions  Moleculaires 
Universite  Pierre  et  Marie  Curie 

N.Gresh 

Laboratoire  de  Bicchimie  Theorique,  associe  au  CMRS 

C,Gary-Bobo,M.Herv4 

Service  de  Biophysique,  Departement  de  Biologie,  CHM  Saclay 


II  existe  de  nombreuses  etudes  experimentales  de 
l'activite  des  polyenes  ionophores  et  de  leur  selectivity 
antif ongique .  L'activite  ionophorique  des  polyenes  est  due 
a  la  formation  de  complexes  avec  les  sterols.  La  stability 
des  complexes  polyenes-sterols  est  fondee  sur  1'energie  de 
deux  types  d' interaction  et  de  liaison: 

1.  les  liaisons  hydrogene  entre  groupements  polaires  du 
polyene  et  1'hydroxyle  sterolique  et  eventuellement 
quelques  molecules  d'eau. 

2.  des  liaisons  Van  der  Waals  entre  la  partie  heptaenique 
du  polyene  et  le  squelette  stero'idique . 

Selon  la  structure  de  la  tete  polaire  du  polyene, 
1'energie  de  la  liaison  hydrogene  est  tres  variable:  si  elle 
est  forte,  le  complexe  ionophorique  est  stable  quel  que  soit 
1 e  sterol;  si  elie  est  faible,  les  interactions  de  Van  der 
Waals  deviennent  predomi  naur.es  e'  il  y  a  alors  selectivite 

v  i  it>  a  /is  Litj  s  oi  tidretV-3  oiei  yio  . 

Une  etude  preaLable  de  la  structure  de  la  tete  polaire  a 
1 ' etat  isole,  ou  en  presence  de  molecules  d'eau,  est  indis¬ 
pensable  avant  1' etude  des  differents  complexes.  Nous  avons 
effectue  des  calculs  par  la  methode  SIBFA:  selon  ce  procede, 
la  "macromolecule"  est  consideree  comme  un  ensemble  de  frag¬ 
ments  relies  par  de  simples  liaisons  et  la  variation 
d'energie  intramoleculai re  est  le  resultat  d'une  somme 
d’energies  d ' interaction  entre  les  fragments.  La  distribution 
de  charges  de  chaque  fragment  est  obtenue  a  partir  d'une 
decomposition  multipolaire  d’une  fonction  d'onde  ab  initio. 
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Modelization  of  the  disulfure  bridge  in 
Proteins  : 

ab-initio- Cl  studies  of  S2H2  and 

Michel  LOOS 

Laborntoire  de  Chi  w  if  Theorique,  U-A.  oLO  C.^RS, 
Uni  versa e  de  Nancy- L  BP  2df). 

F-54506  Vandceuvre-les- Nancy,  Cedex,  France. 


The  conformation  of  the  disulfide  bond  in  proteins  plays  an  important  role  in 
their  properties[lj.  Much  theoretical  work  has  been  done  concerning  the  rotationnel 
barriers  about  the  S-S  bond.  In  this  work,  we  investigated  not  only  the  fundamental 
state  potential  surface  but  also  that  of  the  first  excited  state,  in  order  to  investigate 
the  possible  radiation  d'  mmage  to  the  proteins. 

We  used  two  models  lor  ihe  disulfide  bond  the  SiHi  and  S;j(CHi)2  systems  and 
we  will  report  the  potential  surfaces  around  the  S-S  bound  for  both  of  them  in  the 
fundamental  and  first,  excited  state. 

All  calculations  were  done  with  the  HONDO  progTamm  [21 ,  using  a  4-31G*  basis 
at  the  Cl  level. 

The  disulfide  bonds  are  known  to  exhibit  strong  circular  dichroi'sm.  We  will  also 
report  calculations  of  the  rotatory  strengths  and  dipole  strengths  of  these  systems 
in  the  virtual  orbital  approximation.  We  demonstrate  the  possibility  of  using  inner- 
shell  electronic  spectroscopy  in  the  soft  X-ray  range  to  gain  access  to  the  absolute 
conformation  around  the  disulfide  bond  in  proteins  from  the  circular  dichroi'sm. 
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Analysing  and  Modeling  the  Deformation  of 
Protein  Secondary  Structures 


Catherine  Etchefcest 
Institut  de  Biologie  Physico-Chimique 
13,  rue  Pierre  et  Marie  Curie 
Paris  75005,  France 


Secondary  structures,  most  notably,  alpha  helices  3r.d  beta 
sheets  form  a  considerable  part  of  most  protein  conformations . 
However,  these  structural  elements  rarely  approach  perfect 
helicoidal  symmetry  and  their  participation  in  the  overall 
conformation  often  involves  considerable  deformations.  In  order 
to  attempt  to  understand  the  energetic  nature  of  these  effects  we 
have  made  two  types  of  investigation. 

Firstly,  using  a  new  algorithm  we  have  analysed  the  nature 
of  the  deformations  occurring  in  well  resolved  protein 
structures,  notably  the  bending  of  alpha  helices  and  the  bending 
or  torsion  of  beta  sheets.  Secondly,  we  have  made  theoretical 
calculations  of  deformation  energy  using  a  methodology  which 
enables  us  to  control  easily  the  overall  conformation  of 
oligopeptide  fragments.  We  attempt  to  correlate  the  findings  with 
the  influence  of  the  peptide  sequence  and  to  clarify,  in 
particular,  the  role  of  special  residues  such  as  proline. 
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Describing  Protein  Conformation:  A  New  Mathematical  Approach 


Krystyna  Zakrzewska 
Institut  de  Biolcgie  Physico-Chimique 
13 ,  rue  Pierre  et  Marie  Curie 
Paris  7500d,  France 


Although  a  considerable  number  of  protein  structures  have 
now  been  crystallographically  determined  to  high  resolution  there 
are  still  difficulties  in  extracting  all  the  conformational  data 
that  these  results  contain.  In  particular,  there  is  no  rigorous 
procedure  for  the  precise  location  of  secondary  structures  and 
the  description  of  their  deformations,  or  for  determining  the 
exact  pathway  followed  by  a  folded  polypeptide  backbone. 

An  algorithm,  P-Curves,  recently  developed  in  our 
laboratory  enables  these  goals  to  be  achieved  in  a  clear  and 
general  way.  This  approach  corresponds  to  a  generalized 
helicoidal  description  of  protein  structure  and  yields  both  an 
overall  axis  describing  the  folding  of  the  protein  backbone  and  a 
full  set  of  helicoidal  parameters  locating  each  peptide  in  space. 
The  analysis  of  a  number  of  well  resolved  protein  conformations 
subsequently  allows  us  to  define  the  extent  of  regions  of 
secondary  structure  without  reference  to  local,  and  often 
subjective  measures,  such  as  hydrogen  bonding  or  phi/ps i  angles. 
This  data  also  allows  easy  comparison  of  homologous  protein 
conformations  and  should  facilitate  the  definition  of  sequence- 
structure  correlations. 
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Tautomerism  of  the  Nucleic  Acios  Sases  Revisited: 

From  Isolated  CNon-inter acting)  Bases  to  their  Completes 

Jdsef  S.  Kwi  atkowski ,  Institute  of  Physios.  U.  CcfsoniOMS 

\j&rs  £  t  y »  S“7  ~  i  CC  ror^r. ,  i  ay.d 

Results  are  presented  -from  the  recent  ob  initio 
quantum-mechanical  studies  of  the  relative  stabilities  o-f  the 
isolated  ( nan— i nter acti ng  with  the  molecules  a-f  an  environment) 
tautomers  o-f  the  nucleic  acids  bases.  Calculations  predict  that 
the  "rare"  tautomers  (important  -for  tautomeric  mechanism  of  point 
mutations)  are  higher  in  energy  than  the  "normal"  tautomers  by 
60-65  kJ/mol  -for  adenine  and  thymine  (uracil),  but  only  by  7-3 
kJ/mol  -for  cytosine,  while  in  the  case  o-f  guanine  its  rare 
tautomer  is  lower  by  3  kJ/mol  than  its  normal  tautomer.  The 
predictions  correlate  well  with  recent  experimental  findings  from 
the  infrared  spectroscopic  studies  of  the  bases  isolated  in  inert 
gas  matrices. 

Results  are  also  presented  from  the  combined  approach  (ofc 
initio  calculations  for  relative  stabilities  of  the  isolated 
tautomers  together  with  a  calculation  using  Cl  aver ie’s  method  for 
predicting  i ntermol oecul ar  interaction  energies)  applied  for  study 
the  energetics  of  dimers  and  complexes  consisting  of  different 
tautomers  of  cytosine  and  guanine.  It  is  shown  that  the 
interactions  between  the  bases  change  the  relative  stability  of 
the  individual  tautomers.  particularly  "discriminating"  against 
the  dimers  (complexes)  containing  the  rare  tautomers.  These 
i ntermol ecular  interactions  may  be  an  important  factor  involved  in 
the  complicated  control  (repairing)  systems  of  living  organisms, 
which  discriminates  against  the  rare  (non-canonic)  tautomers  of 
the  bases  during  several  biological  processes. 
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ELECTROSTATIC  BASIS  FOR  PROTEIN  FUNCTIONS:  INCORPORATING  THE  EN¬ 
VIRONMENT  IN  MICROSCOPIC  CALCULATIONS  OF  FREE-ENERGY  AND 
DYNAMICS  OF  CHEMICAL  PROCESSES  IN  SOLUTIONS  AND  PROTEINS 

ARIEH  WARSHEL.  DEPARTMENT  OF  CHEMISTRY,  UNIVERSITY  OF  SOUTHERN 
CALIFORNIA,  LOS  ANGELES,  CA  90089- 1062 

It  has  been  arg  ed  repeatedly  that  electrostatic  free  energies  are  the  key  factors  that 
correlate  structure  ane  function  in  proteins. 1  However,  the  examination  of  this  hypothesis  is 
somewhat  meaningless  without  the  ability  to  perform  microscopic  calculations  of  enzymatic 
reactions  and  related  processes.  Our  developments  in  this  field  during  the  last  15  years  will  be 
reviewed.  This  will  cover  solvent  models  ranging  from  the  simplified,  yet  microscopic, 
Langevin  Dipoles  moder1-2'  to  surface  constraint  all  atoms  models  which  were  implemented 
with  free  energy  perturbation  calculations.^1 3®  The  incorporation  of  microscopic  solvent 
models  in  quantum  chemical  calculations  will  be  considered,  discussing  the  relative 
advantages  of  MO-SCF  approaches®  and  Valence  Bond  approaches.®4 5 6 7^  The  implementation 
of  such  methods  in  calculations  of  activation  free  energies  in  solutions  and  proteins  will  be 
described.  ‘  The  role  of  the  electrostatic  fluctuations  of  the  solvent  in  controlling  dynamical 
effects  will  also  be  discussed  and  the  evaluation  of  dynamical  effects  in  quantum  chemical 
calculations  of  charge  transfer  reactions  will  be  outlined. 
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MONTE-CARLO  FREE  ENERGY  CALCULATIONS  IN  CONFORMATIONAL 
STATISTICS  OF  POLYPEPTIDE  CHAINS 


S.  PREHILAT 

Laboratoire  de  Biophysique  Mol£culaire  -  University  de  Nancy  I 
54506  Vandoeuvre  les  Nancy 


Entropy  and  more  precisely  free  energy  are  the  relevant 
quantities,  besides  the  commonly  evaluated  conformational 
energy,  for  the  theoretical  study  of  macr omolecular  systems 
especially  when  they  are  subject  to  solvent  effects.  But  there 
is  generally  no  direct  access  to  3uch  quantities  because  of 
convergence  problems  in  conventional  Monte-Carlo  simulations. 
Moreover,  even  if  it  is  clear  that  short-range  interactions 
restrict  chain  units  to  a  limited  number  of  low  energy 
conformations,  one  has  to  considere  actually  excluded  volume, 
long-range  electrostatic  interactions  and  solvent  effects  as 
main  factors  for  analysing  the  folding  of  molecular  chains 
into  stable  conformations.  However,  the  introduction  of  such 
long-range  effects  in  calculations  is  a  source  of  major 
difficulties  as  one  can  no  more  express  the  conformational 
energy  as  a  sum  of  terms  associated  to  chain  units  or  first- 
neighbor  pairs. 

However,  the  partition  of  the  total  conformational  energy 
of  a  chain  molecule  into  short-range  and  long-range 
interactions  allowed  us  to  propose  an  effective  new 
approach! 1 )  for  using  Monte-Carlo  calculations  in  statistical 
mechanics  of  molecular  chains.  The  calculation  method  in 
which  the  part  of  energy  coming  from  short-range 
interactions  is  not  recalculated  at  each  step  of  the 
procedure,  is  based  on  the  use  of  statistical  weight  matrices 
for  short-range  terms  and  importance  sampling  for  the  part  of 
energy  associated  to  long-range  interactions. 

Such  an  approach  gives  also  a  way  of  defining  all  the 
statistical  properties  for  reference  states  of  the  molecular 
system.  Therefore,  it  becomes  possible  to  study  the  behavior 
of  the  molecule  under  more  complex  condi t ions ( 2 > .  The  free 
energy  and  entropy  variations  due  to  long-range  interactions 
are  then  calculated.  Moreover,  solvent  effects  can  be 
introduced  by  simple  modifications  of  long-range  terms  of  the 
conformational  energy  as  it  is  generally  done  for  a  potential 
of  mean  force.  Beside  the  determination  of  thermodynamic 
properties  of  the  chain  molecule,  probabilities  of  chain  units 
being  in  their  different  conformational  states  are  also 
estimated.  This  gives  finer  information  on  the  conformations 
of  chain  molecule.  Application  of  these  calculation  methods 
are  made  on  molecular  models  of  peptidic  hormones. 


1.  S.  Pr^milat, J.  Chem.  Phys.  88,3385(1988). 

2.  S.  Pr£milat,J.  Chem.  Phys.  90,5128(1989). 


Dielectric  screening  in  proteins 


Thomas  Simonson1,  David  Perahia2  and  Gerard  Bricogns1 
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proteins  'vnh  "he a-  substrates. 


The  Local  DOianzabiiitv  snourd  van?  tfirouo'hout a  eroiein  nr  men  ■>  ■v.v?  ir  -o  contribute  to- 
its  acuvuv.  To  test  this  hypothesis  -ve  vse  •«  render!  microscopic  model  »c  describe  die  proem's 
dielectric  properties .  end  calculate  'hecreuceiiy  ;ts  susceptibility  in  response  to  e.  set  of  ex*ernel  point 
charges .  The  microscopic  model  desenbes  the  electronic  degrees  of  freedom  by  a  set  of  atomic  point 
polarizabilities  and  the  atomic  positional  degrees  of  freedom  by  normal  mode  dynamics. 


Computer  programs  implementing  me  model  have  been 
the  Charmm  molecular  mechanics  package  (Brooks  at  al,  1983). 
earned  out  on  several  systems  (model  alpha  helices,  cytochrome 


ievelopped  and  incorporated  into 
Numerical  calculations  have  been 
c .  the  disk  of  protein  of  tobacco 


mosaic  '-iris’1  vhich  -end  ro  confirm  our  initial  hypcchesis  ( Simcinson  e*  ai.  1989  ' 
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Calculation  of  the  Free  Energy  of  Intercalation 


by 

B.M.  Pettitt 

Graphite  intercalation  compounds  (or  the  synthetic  metals)  have  many  uses  in  materials  sci¬ 
ence  and  catalysis.  A  statistical  mechanical  theory  is  developed  and  applied  to  study  the 
structural  effects  that  the  thermodynamic  state  of  alkali  ions  have  on  graphite  intercalation 
compounds.  The  system  considered  is  that  of  second  stage  Rb-graphite.  Two  dimensional 
diffraction  patterns  are  computed  and  compared  with  experimental  measurments.  Sensitiv¬ 
ity  to  model  parameters  are  considered.  A  low  order  density  functional  expansion  is  found 
to  adequately  describe  the  interionic  structure  and  free  energies  of  the  system  modeled  as  a 
two  dimensional  one  component  plasma  in  an  anisotropic  external  field.  Work  in  progress  on 
extensions  to  molecular  cases  and  the  calculation  of  interlayer  free  energies  aree  also  presented. 
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Artificial  Intelligence  in  the  Molecular  World 

Marie-Christine  HATON 
CRIN/INRIA 


University  of  Nancy  I 
BP  239 

54506  Vandoeuvre-ies-Nancy  (France) 


Abstract  : 

While  solving  problems  is  a  long  tradition  in  experimental  sciences, 
research  on  how  people  actually  solve  problems  is  relatively  recent.  It  was 
supposed  twenty  years  ago  that  people  knew  how  to  perform  problem 
solving  and  overall  knew  how  to  teach  their  methods. 

This  subject  has  been  recently  deepened  in  the  frame  of  cognitive  sciences, 
namely  by  knowledge  psychologists  and  now  by  knowledge  engineers,  inside 
"Artificial  Intelligence"  (Al).  Al  is  a  branch  of  advanced  information 
processing  dealing  with  these  activities  which  characterize  human 
behaviour  :  knowledge  acquisition  and  structuration,  reasoning,  perception, 
decision  making,  etc.  Al  methods,  tending  to  formalize  human  knowledge 
and  reasoning  processes,  give  us  new  tools  for  a  better  understanding  of 
such  processes. 

Nowadays,  while  the  use  of  computers  is  increasing,  the  amount  and  the 
complexity  of  the  data  involved  in  the  different  domains  of  chemistry 
make  it  imperative  to  discriminate  which  information  is  relevant  and  to 
consider  new  ways  of  processing  it. 

Various  aspects  of  Al  may  concern  the  molecular  world  :  knowledge  and 
reasoning  representation,  heuristic  problem  solving,  planning,  image  and 
natural  language  understanding,  knowledge-based  and  expert  systems, 
"intelligent"  computer-aided  education,  symbolic  learning.  These  aspects 
will  be  discussed,  focusing  on  the  specific  problems  the  chemist  or 
physicochemist  has  to  face  with  (synthesis,  experimental  planning, 
structure  elucidation,  spectra  interpretation,  classification,  information 
retrieval,  fault  diagnosis  in  analytical  instrumentation,  student  education 
and  so  on). 
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IONIC  EFFECTS  ON  BIOflOLECULAR  CONFORMATION 
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MOLECULAR  MODELLING  IN  BIOLOGY, 
with  or  without  QUANTUM  CHEMISTRY. 


Bernard  PULLMAN 


Institut  de  Biologie  Physico-Chimique 
Fondation  Edmond  de  Rothschild 
13,  rue  Pierre  et  Marie  Curie 
75005  Paris. 


Many  people  seem  to  ignore  to  what  an  extent  the 
consideration  or  the  omission  of  certain  fundamental  quantum- 
chemical  concepts  has  played  a  decisive  role  in  the  success  of 
some  and  the  failure  of  others  in  achieving  basic  discoveries 
which  laid  the  foundations  of  molecular  biology. 

Today  when  modelling  procedures  strive  towards  the 
elucidation  of  complex  structures  and  interaction  mechanisms 
the  necessity  of  continuing  to  use  quantum  mechanical  concepts 
and  methods  is  even  greater  than  before.  The  insufficiency  of 
molecular  graphics  based  solely  on  qualitative  considerations 
of  geometrical  fitting,  hydrogen  bonding  capabilities,  and 
similar  pictorial  representations,  will  be  illustrated  in  the 
topical  efforts  of  the  search  for  antituor  DNA  binding  ligands 
with  selective  specificity  for  the  minor  groove  of  AT  or  GC 
base-pair  sequences. 
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CAMBRIDGE. Massachusetts  02 1  38 


HAAS  Verena  Mile 
Univ.  GH  Duisburg 
FB6 

Theoretische  Chemie 

Lotharstrasse  1 

4[00  DUISBURG  1  (RFA) 


KASSAB  E. 

D.I.M. 

Tour  22/ler  etage 
Lmv.  P.  et  M.  Curie 
4,  place  Oussieu 
75252  PARIS  Cedex 


K.OCHANSKI  E.  Mme 
ER  139 

Chtmie  theorique 
Inst,  de  Chimie 

1.  rue  Blaise  Pascal 
B.P.  296/RS 

6700S  STRASBOURG 

KOGELKAJ. 

Lniv.  Rene  Descartes 

Chimie  et  Biochirme  Pharmacoiogiques 

itj,  rue  des  Saints  Peres 

75270  PARIS  Cedex  06 

KWIATKOWSKI  J.S. 

N.  Copernicus  Umv. 

Inst,  of  Phvsics 
S7-100  TORL'N  (Pologne) 

LACROIS-DORE  M.D.Mme 
IBM  FRANCE 
67,  quai  de  la  Rapee 
75012  PARIS 

LAGANT  P. 

Fac.  Pharmacie  Lille 
INSERM  L  279 

3.  rue  du  Professeur  Laguesse 
59045  LILLE  Cedex 

LAHANA  M. 

Aquitaine  Systemes 
Tour  Elf 

2,  place  de  la  Coupole 
9207S  PARIS  LA  DEFENSE 

LAMBERT  C.  Mme 
I.C.I.  PHARMA 
Centre  de  Recherches 
Chemin  de  Willy 
E.l.  la  Pompelle 
B.P.  401 

51064  REIMS  Cedex 

LANGLET  3.  Mme 
Dynamique  des  Int.  Molecuiaires 
Tour  22 
Couloir  22-23 

4,  place  Jussieu 
75005  PARIS 

LANGLET  G. 

CEA/IRDI/DLPC/SCR 

CEN  SACLAY 

91 1 9  1  GIF  SLR  YVETTE 


LAUPRETRE  F.  Mme 

Physico-chimie  structurale  et  macromoleculaire 
ESPCi 

10,  rue  Vauquelin 
75231  PARIS  Cedex  05 

LAVERY  R. 

I.B.P.C 

13.  rue  P.  et  M.  Curie 
75005  PARIS 

LECLERC  3.P. 

Centre  de  recherches  de  Vitrv 
B.P.  1“ 

13,  quai  Jules  Guesde 
94403  VITRY  Si  Seine  Cedex 

LECLERCQ  J.M. 

Lab.  dvnamique  lnt.  molecuiaires 

UPR  AG271  CNRS 

Lniv.P.  et  M.  Curie/Tour  22 

4,  place  Jussieu 

75252  PARIS  Cedex  05 

LEGO AS  M. 

D.C.S.O. 

Ecole  Poly  technique 
91  12S  PALAISEAL 

LEMARECHAL  C. 

1NR1A 
B.P. 105 

7SI53  LE  CHESNAY 

LLUCH  J. 

Depart.  Quimica 

Univ.  Autonoma  de  Barcelona 

OS  193  BELLATERRA  (Barcelona) 

Espagne 

LOOS  M. 

Lab.  Chimie  theonque 
Lniv.  Nancy  1 
B.P.  239 

54506  VANDOELVRE  LE5  NANCY 

MAISSIAT  C. 

Societc  3iomodete 
CEI 

15,  rue  Carnot 
BP  175 

86004  POITIERS  Cedex 
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MAJOUBE  .Vi. 

CEA-DLPC 

CEN  SACLAY 

91  [91  GIF  SLR  YVETTE 

MALLIAVIN  T-.MUe 
DCSO 

Ecole  Polvtechmque 
9112S  PALA1SEAU'  Cedex 


MARK  F. 

Max  Planck  Inst,  fur  Strahlenchemie 
Stiftstr  34-36 
D  a33Q  MULHEIM 
RFA 

MAROLN  R. 

Inc.  Biologique 

Pharmacocairme  moleculaire 
a  av.  de  1’Gbservaroire 
75006  PARIS 

MASSEY  K. 

FPS  COMPUTING 
21.  rue  des  Cevenr.es 
Silic  523 
94633  RUNGI5 

MATHiEX  J.P. 

CRAY  RESEARCH  FRANCE 
7.  rue  de  Tiisitt  75017  PARIS 

.MATHIS  H. 

Chirrne  theorique 
Univ.  Nancy  I 
B.P.  239 

54506  VANDOEUVRE  LES  NANO 

MATT  ALIA  J.M. 

Chimie  inorganique 

Fac.  des  Sciences  St  Jerome 

13397  MARSEILLE  Cedex  13 

MALREL  J.L. 

Procida/Groupe  Roussel-Uclaf 

ST  MARCEL 

130 11  MARSEILLE 

MAZEAU  K. 

CERMAV.CNRS 
B.P.  53x 

38041  GRENOBLE  Cedex 


V1EHANI  S. 

BP  RESEARCH  Center 
Chertsev  Road 
Sunbury  on  Thames 
MIDDLESEX.  LDN.  TW137  LN/UK 

MOY  M. 

AQUITAINE  SYSTEME5 
Tour  Elf 

2.  place  de  la  Coupole 
92078  PARIS  LE  DEFENSE 

MIKOU  A.  Mile 
ICSN/Lab.  R.ViN 
rue  de  ia  Terrasse 
91190  GIF  SUR  YVETTE 

MILLOT  C. 

Lab.  Chimie  theorique 
Univ.  Nancv  1 
B.P.  239 

54506  VANDOEUVRE  LES  NANCV 

MILON  A. 

URA  3i-CNRS 
CENTRE  DE  NEUROCHIMIE 
5,  rue  Blaise  Pascal 
67084  STRASBOURG 

MONTE1L  P. 

CONVEX 

Parc  d'Activites  du  Pas  du  Lac 
Immeuble  le  Daguerre 
9.  av.  Ampere 

78180  MONTICNY  LE  BRETONNEUX 

MONQUE  R. 

INTEVEP  SA 

Apdo  76343  Caracas  1070  A 
Venezuela 

MOREAU  G. 

ROUSSEL  LCLAF 
102.  route  de  Noisy 
93230  ROMAINVILLE 

MORGANTINI  P.Y. 

Depart,  chimie  physique 
Univ.  de  Geneve 
30,  quai  E.  Ansermet 
1211  GENEVE  4  (Suisse) 

MOR1ZE  I.  Mme 
Rhone  Poulenc  Sante 
Modelisation  moleculaire 
1  3.  rue  Jules  Guesde 
94403  V1TRY  S.  Seine  Cedex 


MULLER  N.  Mme 
Inst.  Henri  Beaufour(IHB) 

17,  av.  Descartes 

92350  LE  PLESS1S  ROBINSON 

NICKLAS  K. 

Inst,  fur  Phvsikalische  Chemie 
TH  DARMSTADT 
Petersenstr.  20 
6100  DARMSTADT  (RFA) 

OHLENBUSCH  H. 

Inst,  ae  Physique  biologique 

4,  rue  Kirschleger 

670S5  STRASBOURG  Cedex 

PANTEL  G. 

Rhone  Poulenc  Same 
Centre  de  recherches  de  Vitrv 
BP  14 

13.  quai  Oules  guesde 

944Q3  VITRY  SLR  SEINE  Cedex 

PEPE  G. 

CRMC2-CNRS 
CAMPUS  LUMINY 
Case  9 1 3 

132SS  MARSEILLE  Cedex  9 

PESQUER  M. 

Phvsicochimie  theorique 
CNRS-UA  503 
Lniv.  Bordeaux  I 
33405  TALENCE  Cedex 

PETTITT  B.M. 

Chemistry  department 
Univ.  of  Houston. 

4800  CALHOUN 

HOUSTON  TEXAS  77204-5641 

USA 

PIRIOL  3.M. 
Biochimie-Enzvmologie 
Inst.  Gustave  Roussy 
39,  rue  Camille  Desmoulins 
94805  VILLE0UIF  Cedex 

POLTEV  V. 

Inst,  of  Biological  Physics 
USSR  Academy  of  Sciences 
Pushchmo,  Moscow  Region 
142292  USSR 


PONCIN  M. 

IBPC 

13,  rue  P.  et  M.  Curie 
75005  PARIS 

POTHIER  J. 

Biologie  Enzvmologie 
Inst.  Gustave  Roussy 
39.  rue  Camille  Desmoulins 
94805  VILLE3UIF  Cedex 

PRE.XiILAT  S. 

Biophysique  molecuiaire 
Lniv.  Nancv  1 

54506  VANDOEUVRE  LES  NANCY 

PULLMAN  A. Mme 

Inst,  de  Biologie  Phvsico  Chimique 

13.  rue  P.  et  M.  Curie 

75005  PARIS 

PULLMAN  B. 

Inst,  de  Biologie  Phvsico-Chimique 
13.  rue  P.  et  M.  Curie 
7  5005  PARIS 

PLOCKY.N  M. 

Jobin  Yvon 

16-18  rue  du  Canal 

91160  LONGJUMEAU 

RINALDI  D. 

Chimie  theorique 
Univ.  Nancv  I 
B.P.  239 

54506  VANDOEUVRE  LES  NANCY 

R1VAIL  3.L. 

Lab,  Chimie  theorique 
Unit.  Nancv  ) 

B.P.  239 

54506  VANDOEUVRE  LES  NANCY 

ROBERT  B. 

CENS/SBPH 

CEN  SACLAY 

91  191  GIF  SLR  'r  VETTE 

RUBIN  C.Mlle 
Inst.  Chimie 
1,  rue  Blaise  Pascal 
6700  STRASBOURG 
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RUIZ  M. 

Chimie  theorique 
Univ.  Nancv  I 
B.P.  239 

59506  VANDOEUVRE  LES  NANCY 

RULLMANN  J.A.C. 

Department  of  Organic  Chemistry 

Univ.  of  Utrecht 

Paduataan  S 

3 5 S A  CH  UTRECHT 

(Pays  Bas) 

RAHMOLN1  A. 

Chimie  theorique 
Inst,  de  Chimie 
1.  rue  Biaise  Pascal 
6700  STRASBOURG 

REPELIN  Y.  \lme 
ECOLE  CENTRALE 
Grande  vuie  des  vignes 
92290  CHATENAY  MALABRY 

SAARMETS 
Synthelabo  Recherche 
31,  av.  P.  Vaillant  Couturier 
92200  BAGNEUX 


SOLMAJER  T. 

Boris  Kidric  Chemistry  Inst. 
Hajdrihova  19 

6100  LJUBLJANA  SLOVENIA 
Y  ougoslavie 

SOULIE  E. 

IRDI/DESICP/DIPC/SCM 

CEN  SACLAY 

91191  GIF  SUR  YVETTE 

SCUMPAS15  D.M. 

Max  Planck  Inst  fut  Bioph.  Cnemistry 
p.o  Box  2891 
3900  GOTTINGEN  (RFA) 
et  LOS  ALAMOS  nat.  lab. 

LOS  ALAMOS  $7595 
USA 

ST  AW  ARZ  B. 

Centre  de  Biophvsicue  Moleculaire 
1A.  av.  de  la  Recherche  Scientitiaue 
95071  ORLEANS  Cedex  2 

STONE  A. J. 

Univ.  Chemical  Lao. 

Lensfietd  Road 

CAMBRIDGE  CB2.  1  EW  (CB) 


SCHLENKR1CH  M. 

Inst,  fur  Physikalische  Chemie 
TH  Darmstadt  BRD 
Pecersenstr  20 
6100  DARMSTADT.  BRD 

SCHVVAAB  F. 

CRIN 

Univ.  Nancv  1 

59506  VANDOEUVRE  LES  NANCY 

SHUSTOROVICH  E. 

CORPORATE  RESEARCH  LAB 
Eastman  Kodak  Company 
CORL  B-S 1 

ROCHESTER. NY  1 9650-0200  HU.S.A.) 
SIRI  D. 

CRMC2-CNRS 
Campus  Luminy 
Case  913 

13288  MARSEILLE  Cedex  9 


SMITH  J. 

CISI-CEA 

Service  de  Biophysique 
Depart,  de  Biologie 
9119  1  GIF  SUR  YVETTE 


SUN  Jian-Sheng 
Biophvsique.lnserm  201 
CNRS  LA  981 

MUSEUM  NATIONAL  D'HISTOIRE  NATURELLE 
93,  rue  Cuvier 
75031  PARIS 

SURCOUF  E.  Mme 
RHONE  POULENC  SANTE 
Centre  de  recherches  de  Vitrv 
BP  |9 

1  3  quai  Jules  Guesde 
99903  VITRY  sur  Seine  Cedex 

TAYLOR  J. 

BiOSY'M  Technologies  Inc. 

15  av.  Victor  Hugo 
751  16  PARIS 

TORRENS  F. 

Departament  de  Quimica  Fisica 
Univ.  ■‘t  Valencia 
Dr.  Moliner  50 
96100  BURJASSOT  (Espagne) 

TRINQU1ER  G. 

Physique  quantique 
Univ.  Paul  Sabatier 
118,  route  de  Narbonne 
31062  TOULOUSE  Cedex 
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TROYANOWSKY  C. 

DIVISION  DE  CHIMIE  PHYSIQUE 
10,  rue  Vauquelin 
75005PARIS 

URBAIN  D. 

FPS  COMPUTING 
21.  rue  des  Cevennes 
Si  Lie  523 
94633  RUNG'S 

V  A  ISM  AN  1. 1. 

Inst.  of  non  Aqueous  Solutions 
Chemistry  of  the  USSR 
Academy  of  Sciences 
Akademicheskaka  Street  1 
153045  IVANOVO 
USSR 

VAN  GU5TEREN  W.F. 

Depart,  of  Physical  Chemistry 
Univ.  of  Groningen 
Nvenborgh  16 
9747AG  GRONINGEN 
Pa>s  Bas 

VERGELATI  C. 

RHONE  POULENC  RECHERCHES 
Centre  de  recherches  des  Carrieres 
S5  av.  des  Freres  Perret 
B.P.  62 

69192  SAINT  FONS  Cedex 

VERGOTEN  G. 

Fac.  de  Pharmacie  Lille 
INSERM  U279 

3.  rue  du  Professeur  Laguesse 
0904  >  LILLE  Cedex 

V1G.NE  F.  Mme 

1. R.C. 

2, av.  Albert  Einstein 

69626  VILLEURBANNE  Cedex 

VOISIN  C. 

Chimie  theorique 
Univ.  de  Nancv  I 

B. P.  239 

54506  VANDOEUVRE  LE5  NANCY 

VOLTZ  R. 

C. R.N. 

Univ.  Louis  Pasteur 
BP  20 

67037  STRASBOURG  Cedex 


WARSHEL  A. 

U.  Southern  California,  Chemistry, 

U.  Park, 

LOS  ANGELES, CA  90CS9-1G61 
U.S.A. 

WTPFF  G. 

RMN  et  Modelisation  moiecuiaire 
Inst,  de  Chimie 
1.  rue  B.  Pascal 
6700  STRASBOURG 

ViODAK  S.Mme 
Univ.  Libre  de  Bruxelles 
UCMB 
CP  160 

Av.  Paul  Heger  P2 
B  1050  BRUXELLES 

WURFLINGER  A. 

Lehrstuhl  fur  Physikalische  Chemie  II 
Ruhr  Univ.  Bochum 
Universitatsstrasse  160 
46  30  BOCHUM  1 

ZAKRZEWSKA  K.  Mme 
I.P.B.C. 

1  3,  rue  P.  et  M.  Curie 
7  5005  PARIS 

ANGYAN  3.G. 

Spectrochimie  moiecuiaire 
Univ.  P.  et  \l.  Curie 

4,  place  Jussieu 
75005  PARIS 

BOUFERGUENE  A. 

Chirr, ic  r  rique 

Univ.  .  _v  1 
B.P.  239 

54506  VANDOEUVRE  LES  NANCY 

SAM  AM  A  M. 

Biostructures 

5.  rue  Gustave  Adolphe  Hirn 
67000  STRASBOURG 


